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Abstract 

This  is  a  brief  sunmiary  of  a  study  to  reduce  di^raction  from  an  infmite  two 
dimensional  edge  using  tapered  periodic  surfaces.  The  concept  is  to  use  the  irrq)edance 
characteristics  of  a  periodic  surface  to  create  a  tapered  impedance,  much  like  a  tapered 
resistance  commonly  used  for  di^action  reduction.  This  approach  has  several  advantages 
over  tapered  resistive  treatments  to  include  ease  of  fabrication,  betto-  engineoing  control  of 
the  impedance  and  negligible  loss  absorbed  by  the  taper.  In  this  study,  four  periodic  tapers 
were  design,  built  and  tested.  These  tapers  were  also  modeled  using  a  variety  of  numoical 
techniques.  The  results  from  this  study  were  quite  promising. 

The  tapers  were  designed  using  a  periodic  moment  method  (PMM)  code. 

Assuming  local  periodicity  the  PMM  code  was  used  to  determine  the  impedance  of  doubly- 
infinite  structures.  The  taper  was  developed  by  determining  the  best  taper  response  for  a 
given  bandwidth. 

An  antenna  structure  to  simulate  near  grazing  incidence  of  a  plane  wave  on  a  two 
dimensional  knife  edge  was  used  to  test  the  taper  designs.  The  uq)ers  were  then  built  and 
measured  at  the  compact  range  at  OSU-ESL  over  the  frequency  range  2-18  GHz. 

The  tapers  were  also  modeled  using  several  different  numerical  techniques.  The 
first  designs  were  modeled  using  method  of  moments  (MoM)  techniques.  The  srcond 
capacitive  design  was  modeled  using  the  finite  difference  time  domain  (FDTD)  technique. 

The  results  from  the  measurements  were  positive.  The  results  from  the  numerical 
calculations  supported  the  measurements  and  displayed  the  bandwidths  indicated  from 
earlier  PMM  calculations.  The  inductive  designs  developed  were  not  very  successful ,  but 
the  capacitive  taper  proved  to  be  quite  effective  over  a  wide  frequency  range. 
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A  STUDY  OF  TWO  DIMENSIONAL  TAPERED  PERIODIC  EDGE 
TREATMENTS  TO  REDUCE  WIDEBAND  EDGE  DIFFRACTION 


Chapter  I:  Objectives  of  the  Study 


Introduction 

Overview.  The  purpose  of  this  chapter  is  to  describe  a  study  of  an  edge 
diffraction  reduction  technique  using  a  tapered  periodic  geometry.  Many  efforts  have 
previously  been  made  to  reduce  edge  diffraction;  however,  these  techniques  did  not  use 
tapered  periodic  surfaces.  Specifically,  the  purpose  of  this  study  is  to  examine  the  use  of 
tapered  periodic  structures  to  reduce  wideband  edge  diffraction  near  grazing  incidence. 
These  tapers  will  have  specific  application  to  antenna  design  for  high  power,  time  domain, 
and  wideband  antennas. 

The  first  section  of  this  chapter  describes  diffraction  and  the  undesirable  effects  of 
diffraction  in  antenna  design.  The  second  section  describes  the  general  problem  statement 
and  explains  the  concepts  proposed  for  reducing  diffraction  for  both  principle  polarizations. 
The  third  section  describes  the  current  technologies  to  reduce  diffiraction.  The  last  section 
describes  the  assumptions  made  in  the  study. 

Background 

Purpose.  The  purpose  of  this  section  is  to  describe  diffraction  and  the  undesirable 
effects  of  diffraction  in  antenna  design.  This  section  defines  terms  associated  with 
diffraction  and  antenna  design  used  throughout  this  paper. 
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Diffractidn.  Diffraction  is  an  important  topic  of  study  in  electromagnetism.  If  an 
incident  field  illuminates  a  scatterer,  then  the  total  electric  field  in  a  region  is  defined  as  the 
incident  plus  the  scattered  fields  (diffracted  and  reflected  fields).  In  the  far  field  region,  the 
reflected  field  is  present  only  at  discrete  angles  for  a  flat  scatterer  (see  Figure  1 .  Region  2). 
Diffraction,  however,  contributes  scattered  fields  at  all  angles  (see  Figure  1,  Regions  1-3), 
which  makes  diffraction  the  most  prominent  scattering  mechanism  for  far  field 
measurements.  The  shadowed  region  (see  Figure  1,  Region  3)  fields  consist  of  diffracted 
fields  only.  Edge  diffraction  is  defined  as  diffraction  coming  from  a  sharp  edge  of  the 
scatterer.  Knife  edges  are  edges  which  are  made  by  an  infinitesimally  thin  plates.  This 


work  will  concern  the  study  of  reducing  knife  edge  diffraction. 
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Figure  1.  Diffraction  from  a  Conducting  Strip.  This  is  a  2D  conducting 
strip  illuminated  by  an  incident  plane  wave.  The  two  scattering  mechanisms  are 
reflection  and  diffraction. 
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Antenna  Design.  A  major  use  of  the  study  of  diffraction  is  to  optimize  antenna 
design.  Several  parameters  determine  the  quality  of  an  antenna.  Some  of  these  factors  are 
directivity,  gain,  ringing,  and  pulse  distortion  [1].  Diffraction  from  the  edges  of  an  antenna 
adversely  effects  all  of  these  factors.  Therefore,  reduction  of  diffraction  from  the  edge  of 
antennas  has  been  a  main  source  of  study  in  electromagnetism.  Directivity  is  a  measure  of 
how  much  of  the  energy  of  the  signal  is  contained  in  the  boresight  view  of  the  antenna. 
Diffraction  causes  sidelobes  and  reduces  the  directivity  of  the  antenna.  Gain  is  the  measure 
of  the  boresight  transmission  of  the  antenna  in  relation  to  the  input  signal  to  the  antenna. 
Diffraction  spreads  the  energy  out  over  all  angles  and  therefore  reduces  the  gain  of  the 
antenna.  Ringing  is  an  undesirable  effect  caused  by  the  diffraction  from  the  edge  of  an 
antenna  radiating  back  into  the  antenna  and  causing  periodic  residual  radiation.  Ringing 
can  give  false  signals  for  both  transmitting  and  receiving  transient  antennas.  Pulse 
distortion  is  any  change  in  the  radiated  wave  from  the  input  pulse.  Pulse  distortion  is 
caused  by  both  geometry  and  materials  having  time  dispersive  and  frequency  dependent 
properties.  Ringing  and  pulse  distortion  are  of  concern  when  working  with  time  domain 
antennas.  Reduction  in  diffraction  from  the  edges  of  the  antennas  will  have  a  positive 
affect  on  all  of  the  mentioned  parameters  (see  Figure  2  for  applications). 

Radar  Cross  Section.  Another  advantage  of  reducing  diffraction  is  to  reduce 
radar  cross  sections  (RCS)  of  targets.  Radar  cross  section  is  a  measure  used  to  quantify 
the  magnitude  of  the  returned  signal  of  a  target  at  given  angles  and  frequencies.  Since 
reflection  only  accounts  for  returns  at  discrete  angles  (assuming  flat  surfaces),  reduction  in 
diffraction  is  of  high  concern.  Presently,  materials  are  used  with  high  magnetic  loss  to 
reduce  diffracted  fields.  These  materials  are  quite  expensive  and  are  generally  effective  on 
only  a  small  frequency  range.  The  tapers  proposed  here  would  also  be  of  use  to  reduce 
traveling  waves,  which  would  be  helpful  when  trying  to  make  scattering  measurements  and 
measuring  returns  only  from  certain  aspects  of  the  structure  (see  Figure  2  for  applications). 
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Trailing  Edge  Diffraction  Reduction 
Figure  2.  Possible  Uses  of  Geometrical  Tapers.  This  diagrams  shows 
possible  uses  of  the  geometrical  tapers  proposed  in  this  study.  The  U  pers  could 
possibly  be  used  to  reduce  trailing  edge  diffraction  or  diffraction  from  the  edge  of 
antennas. 


Problem  Statement 

Purpose.  The  purpose  of  this  section  is  to  define  the  problem  to  be  studied  in  this 
paper.  This  section  will  first  outline  the  objectives  for  understanding  time  domain  edge 
diffraction.  This  section  will  then  detail  the  concept  of  the  periodic  tapers  to  reduce 
diffraction. 

Determine  Configurations  For  Edge  Tapers.  The  first  step  in  this  study  is 
to  determine  what  configurations  create  the  tapered  effect  desired.  The  concept  is  to  create 
an  impedance  taper  by  varying  the  geometry  of  the  periodic  structures  to  impedance  match 
the  ground  plane  to  free  space.  The  work  presented  is  directed  at  creating  periodic  surfaces 
which  have  total  reactive  impedance  and  no  resistive  impedance.  Resistive  impedance 
tapers  absorb  energy  (through  heat  dissipation)  and  would  not  be  useful  for  high  powered 
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microwave  applications.  If  the  impedance  is  kept  either  capacitive  (negative  reactance)  or 
inductive  (positive  reactance),  then  the  phase  shift  will  not  result  in  time  dispersion  (only  a 
time  shift),  because  the  entire  scattered  field  would  be  shifted  by  the  same  phase  shift  (±n). 
The  first  tapers  were  designed  to  not  have  a  phase  shift,  but  to  increase  bandwidth  and 
performance  this  criteria  was  not  used  in  the  second  design.  The  designs  studied  in  this 
paper  were  determined  by  finding  geometries  which  are  completely  reactive  and  could  be 
tapered  over  a  large  frequency  range.  The  designs  were  initially  developed  using  a  Periodic 
Moment  Method  (PMM)  code  [2].  This  code  was  used  to  give  information  about  the 
bandwidth  and  effectiveness  of  the  design. 

Experimentally  Verify  the  Effectiveness  of  the  Tapered  Surfaces.  The 
designs  developed  from  PMM  were  then  tested.  The  designs  were  mounted  on  a  platform 
with  a  long  knife  edge  and  then  the  difft'acted  fields  were  measured  in  the  shadow  region. 

If  the  tapers  are  effective,  there  should  be  a  noticeable  reduction  in  the  shadow  region 
fields.  To  match  the  idea  of  reduction  of  diffraction  in  antenna  design,  the  diffraction  from 
near  grazing  incidence  was  measured.  The  procedures  and  results  from  these 
measurements  are  presented  in  more  detail  in  future  chapters. 

Computationally  Validate  the  Results  After  the  tapers  were  measured,  the 
next  step  was  to  see  if  computational  tools  exist  to  properly  model  this  measurement  For 
the  first  designs.  Method  of  Moments  (MoM)  codes  existed  to  measure  the  reduction  in 
diffraction.  The  MoM  codes  available  could  not  model  the  second  designs,  since  these 
structures  involved  dielectrics.  A  Finite  Difference  Time  Domain  (FDTD)  code  [3]  was 
used  to  validate  the  measurements  of  the  second  capacitive  design.  No  code  was  available 
that  had  accurate  solutions  for  the  second  inductive  taper. 
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Current  Technologies 

Purpose.  This  section  details  some  of  the  current  technologies  used  to  study  and 
reduce  diffraction.  More  detail  on  these  technologies  may  be  found  in  the  next  chapter. 

Resistive  Loading  for  Time  Domain  Antennas.  There  have  been  several 
papers  on  resistively  loading  time  domain  antennas  to  reduce  diffraction  and  thereby 
improve  the  performance  of  the  antenna.  Maloney  et  al  [4]  described  a  method  for 
optimizing  the  resistive  loading  of  a  time  domain  antenna .  Maloney's  antenna  was  a 
biconical  antenna  which  is  analogous  to  a  double  wedge  in  2D.  This  antenna  was 
resistively  loaded  such  that  the  impedance  of  the  source  was  slowly  matched  to  the 
impedance  of  free  space.  This  antenna  design  was  extremely  effective  for  reduction  in 
ringing  and  for  reduction  of  sidelobes  and  edge  diffraction  for  a  transient  (time  domain) 
pulse.  However,  this  antenna  design,  like  all  other  resistively  loaded  antennas,  was  not 
efficient,  over  14%  of  the  energy  was  dissipated  in  the  loading  which  causes  severe 
problems  when  trying  to  use  antennas  for  high  power  applications  such  as  impulse  radar  or 
high  powered  microwave  weapons.  A  resistively  loaded  antenna  would  melt  from  the  heat 
of  the  dissipation  for  these  applications.  Oearly,  methods  other  than  resistive  loading  need 
to  be  developed  for  use  in  high  powered  systems. 

Thin  Wire  Measurements.  Thin  wire  configurations  have  been  studied 
extensively  in  the  frequency  domain.  The  frequency  domain  analysis  has  been  mostly 
concerned  with  thin  strip  application  to  frequency  selective  surfaces  (FSSs).  These 
surfaces  are  useful  for  radomes  and  other  structures  where  only  a  given  frequency  range  is 
of  interest  and  other  frequencies  would  be  a  nuisance  or  even  a  threat  The  study  of  these 
structures  has  principally  focused  on  doubly-infinite  arrays.  The  information  that  is  most 
useful  to  this  study  is  how  FSSs  can  be  truncated  using  edge  treatments  for  a  given 
frequency  range. 
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Thin  wire  configurations  were  also  studied  as  a  singular  subject.  Analytical 
solutions  were  developed  by  Wait  [5]  and  Yoder  [6].  These  techniques  are  quite  simple  to 
implement  and  give  the  scattering  from  an  infinite  array  of  infmitely  long,  thin,  parallel 
wires.  Others  have  studied,  singly  and  doubly  periodic  structures  both  in  free  space  and  in 
the  presence  of  a  dielectric.  Felsen  and  Carin  [71[81[9]  authored  several  papers  where  they 
demonstrated  the  use  of  Floquet  mode  analysis  for  the  study  of  2D  thin  strips.  One 
important  result  of  the  research  conducted  by  Felsen  and  Carin  is  that  they  developed  a 
method  for  applying  Roquet  mode  analysis,  which  is  generally  for  infinite  periodic 
structures,  to  finite  non-periodic  structures.  This  analysis  does  not  lend  itself  very  well  to 
the  study  presented  in  this  paper,  because  this  technique  does  not  handle  dielectrics. 

Periodic  Method  of  Moments.  To  study  FSSs,  researchers  often  use  the 
PMM  code  mentioned  earlier  [2].  This  PMM  code  can  model  doubly  infinite  periodic 
structures.  This  code  can  include  dielectrics,  networks,  and  work  over  a  wide  range  of 
frequencies.  Since  each  element  in  the  doubly  infinite  array  is  identical,  this  PMM  code 
cannot  be  used  to  model  the  periodic  tapers.  An  in-depth  discussion  of  the  workings  of 
PMM  and  the  code  developed  by  Henderson  is  presented  in  the  next  chapter. 

Method  of  Moments  Numerical  Modeling.  One  of  the  most  accurate 
solutions  for  frequency  domain  electromagnetics  is  the  MoM.  MoM  uses  basis  and  testing 
functions  to  determine  the  equivalent  scattering  currents.  These  equivalent  scattering 
currents  can  then  be  transformed  into  the  scattered  fields  using  the  radiation  integral.  The 
MoM  codes  available  can  only  handle  perfect  electrically  conducting  (PEC)  materials  or 
PEC  loaded  materials  of  finite  extent,  not  lossless  dielectrics.  This  numerical  modeling 
technique  was  used  to  validate  some  of  the  impedance  designs,  but  another  technique  was 
necessary  for  the  designs  which  included  dielectrics. 
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Finite  Difference  Time  Domain  Numerical  Modeling.  There  are  several 
* 

techniques  used  to  model  time  domain  electromagnetics.  A  popular  numerical  nKxleling 
approach  is  to  solve  Maxwell's  equations  directly  in  the  time  domain  using  FDTD.  FDTD 
handled  the  complex  structures  involved  in  the  two  dimensional  designs  that  involved 
dielectrics.  The  specifics  of  this  technique  are  given  in  more  detail  in  the  following  chapter. 

Assumptions 

Purpose.  The  purpose  of  this  section  is  to  define  the  major  assumptions  made  in 
this  study.  These  assumptions  are  made  for  clearer  understanding  of  the  results.  Some 
assumptions  have  different  impacts  on  analytical,  numerical  and  experimental  accuracy. 
Each  part  of  this  section  describes  the  assumptions  made  and  the  consequences  of  each 
assumption.  In  Chapter  4,  the  effects  of  these  and  other  assumptions  ntade  will  be 
explained  in  detail  and  contrasted  to  the  results  from  the  study. 

Plane  Wave  Incidence.  The  incident  wave  was  assumed  to  be  a  plane  wave.  A 
plane  wave  is  a  wave  which  has  the  same  field  values  for  every  point  in  a  plane  at  a  given 
time.  Generally,  if  the  incident  source  is  far  from  the  scatterer  (more  than  5  wavelengths), 
the  plane  wave  approximation  is  accurate.  This  assumption  is  useful  because  plane  waves 
are  the  easiest  waves  to  numerically  and  analytically  evaluate.  The  platform  used  in  the 
measurements  was  r.  jfficiently  large  enough  that  the  source  was  several  wavelengths  from 
the  long  knife  edge  (5X  at  the  lowest  frequency). 

PEC  and  Lossless  Dit*pctrics  Materials.  In  the  initial  study,  the  periodic 
structures  were  constructed  completely  out  of  PEC  material.  Beginning  the  stiKly  with 
PEC  material  had  two  major  advantages:  PEC  material  does  not  dissipate  energy  and  PEC 
has  the  same  infinite  conductivity  over  all  frequencies.  The  frequency  independence 
prevented  pulse  distortion  due  to  the  materials  and  made  the  pulse  distortion  be  a  function 
of  geometry  only.  Experimentally,  materials  with  near  infinite  conductivities  (copper. 
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aluminum)  were  readily  available  to  build  the  structures.  For  FDTD  aiKl  MoM  numerical 
analysis,  non-PEC  material  could  be  added  easily  in  future  studies.  Unfortunately,  when 
trying  use  particular  resistive  values  in  a  taper  design,  materials  may  not  be  readily  available 
which  have  the  proper  impedance  value.  This  problem  will  make  experimental  validation 
of  non-PEC  materials  and  resistive  tapers  much  more  difficult  For  the  frequencies  of 
interest  the  dielectrics  used  in  this  study  were  nearly  lossless. 

Structures  Are  Two-Dimensional  (2D).  To  study  the  2D  case,  the  entire 
configuration  was  considered  to  be  infinite  and  constant  in  the  z  direction.  For  both  FDTD 
and  MoM  numerical  techniques,  2D  analysis  was  much  faster  and  required  less  memory 
than  three  dimensional  (3D)  analysis.  Therefore,  the  2D  assumption  allowed  local, 
inexpensive  computational  facilities  to  be  used  for  the  analysis.  Experimentally,  a  2D  wire 
is  impossible  to  realize,  but  if  the  wires  are  made  at  least  5  times  longer  than  a  wavelength 
at  the  lowest  frequency,  the  2D  approximation  is  accurate  within  experimental  error  for 
predicted  patterns  in  the  transverse  plane.  Another  problem  experimentally  is  that  the 
structures  must  be  kept  straight  and  constant  in  the  z  direction.  To  solve  this  problem,  the 
tapers  were  attached  to  styrofoam  for  stability  and  alignment  Styrofoam  is  an  easy 
material  to  work  with  and  has  an  index  of  refraction  near  unity.  The  2D  assumption  is  the 
one  assumption  that  should  present  the  most  problems  when  trying  to  compare  numerical 
data  to  experimental  data.  Further  explanation  of  construction  of  the  2D  configuration  will 
be  presented  in  Chapter  3. 

Loco/  Periodicity.  Since  the  structures  measured  were  almost  periodic,  methods 
of  modeling  for  periodic  structures  were  used  in  the  design  process.  To  use  a  completely 
periodic  method  for  these  tapered  structures,  local  periodicity  must  be  assumed.  Local 
periodicity  means  that  the  changes  in  the  immediate  area  of  each  scatterer  are  so  small  that 
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scaRering  can  be  approximated  using  periodic  techniques.  If  the  tapers  are  large  enough, 
then  the  local  changes  will  be  small  and  the  periodic  techniques  can  be  used  to  nKxlel  the 
tapers  accurately. 

Methodology 

Purpose.  The  purpose  of  this  section  is  to  outline  the  methodology  used  in  this 
study.  Each  subsection  describes  a  different  stage  of  the  study.  Also  included  in  each 
subsection  are  important  assumptions  and  considerations  for  error  that  must  be  taken  into 
account. 

Design  the  Tapers  Using  PMM.  The  Hrst  step  was  to  develop  the  taper 
designs.  To  do  this,  reflection  and  transmission  coefficients  were  calculated  from  periodic 
designs  using  Henderson's  PMM  Code  [2].  The  ends  of  the  taper  designs  were  developed 
by  trying  to  match  the  ground  plane  and  the  free  space  edges  of  the  taper.  The  designs 
were  then  geometrically  tapered  in  a  linear  fashion  from  one  edge  to  another.  The  results 
from  this  part  of  the  study  were  then  transformed  into  equivalent  impedances.  These 
impedance  tapers  are  displayed  as  impedance  versus  length  or  frequency  in  Chapter  4. 

Experimentally  Confirm  the  Results  in  the  Frequency  Domain.  After 
the  designs  were  developed,  the  next  stage  was  to  experimentally  confirm  the  numerical 
results.  A  platform  was  constructed  to  measure  shadow  region  diffraction  from  a  grazing 
incidence  plane  wave.  This  structure  was  measured  both  bare  and  with  the  periodic  tapers 
at  the  Ohio  State  University  Electroscience  Laboratory  (OSU-ESL)  compact  range. 

Because  compact  ranges  use  parabolic  reflectors,  the  far  field  observation  point 
approximation  for  a  compact  range  was  accurate.  The  measurements  of  most  concern  were 
comparisons  of  azimuth  scans  of  the  bare  structure  to  tapered  structures.  The  scans  of  the 
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bare  structure  were  compared  to  the  Uniform  Theory  of  Diffraction  (UTD)  results  for 
validation.  Also  frequency  sweeps  at  discrete  angles  were  taken  to  provide  wideband 
frequency  information. 

Numerically  Validate  the  Results.  The  next  step  was  to  use  numerical 
techniques  to  validate  the  experimental  results.  The  analytical  results  were  validated  using 
two  separate  types  of  codes.  The  first  type  of  code  was  a  MoM  code.  Two  separate  MoM 
codes  were  used  to  model  the  first  impedance  taper  designs.  These  two  codes  handled  the 
two  polarizations ‘used  in  this  study.  The  second  type  of  code  was  a  scattered  field  FDTD 
code  by  Dr.  Ray  Luebbers  [3].  This  code  was  purchased  by  the  U.  S.  Air  Force  and  has 
been  experimentally  validated  for  many  different  geometries.  The  Luebbers'  code  was  run 
in  a  matter  of  a  several  hours  on  a  Silicon  Graphics  machine. 

Summary 

Overview.  This  chapter  has  outlined  the  study  to  be  presented  in  this  thesis.  The 
first  section  gave  a  general  overview  to  the  background  and  reasons  for  reducing 
diffraction.  The  second  section  outlined  the  problem  to  be  studied  in  this  thesis.  The  third 
section  described  current  technologies  used  to  reduce  diffraction.  The  next  two  sections 
gave  the  assumptions  taken  in  this  paper  and  the  scope  of  the  study.  The  last  section  gave  a 
brief  outline  of  the  methodology  to  be  taken  in  researching  this  topic.  More  detail  on  all  of 
these  topics  is  presented  in  the  following  chapters. 
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Chapter  2:  Summary  of  Current  Knowledge 


Introduction 

Overview.  This  chapter  overviews  current  knowledge  that  pertains  to  the  study. 
The  first  section  covers  current  knowledge  of  an  infuiite  parallel  wire  grids.  The  second 
section  covers  meanderline  designs.  The  third  section  explains  the  abilities  and  limitations 
of  numerical  techniques  to  include  the  Uniform  Theory  of  Diffraction  (UTD),  the  Method 
of  Moments  (MoM),  the  Periodic  Moment  Method  (PMM),  and  the  Fmite  Difference  Time 
Domain  (FDTD).  Each  section  covers  how  the  knowledge  obtained  from  each  author  can 
be  used  in  the  design  and  understanding  of  the  tapers. 

Inductive  and  Capacitive  Periodic  Structures 

Overview.  There  are  several  different  types  of  periodic  structures  that  produce 
either  capacitive  or  inductive  impedance.  In  this  study,  the  inductive  properties  of  parallel 
wire  grids  and  meanderlines  are  used  for  the  E-field  orthogonal  to  the  plane  of  incidence. 
The  capacitive  prpperties  of  thin  strips  are  used  for  the  E-field  parallel  to  the  plane  of 
incidence.  To  increa.se  the  continuity  of  the  taper,  the  meanderline  design  was  changed  into 
an  skewed  slot  array  grid  (by  decreasing  the  spacing  between  elements).  Analytical 
solutions  did  exist  for  the  thin  strip  and  thin  wire  grids,  but  numerical  modeling  of  these 
structures  was  far  more  accurate  and  faster  than  analytical  techniques. 

Parallel  Wire  Grids  and  Thin  Strip  Gratings 

Overview.  Parallel  wire  grids  have  been  studied  for  several  years.  Several 
analytical  and  numerical  solutions  for  parallel  wire  grids  have  been  developed.  A  parallel 
wire  grid  was  used  for  the  first  inductive  taper  where  the  E-field  is  polarized  parallel  to  the 
wires.  Due  to  restrictions  in  manufacturing  capability,  a  constant  wire  radius  was 
assumed.  To  change  the  inductance,  the  spacing  between  the  wires  was  changed.  Wires 
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spaced  close  to  each  other  displayed  low  inductance  (ground  plane  side)  and  wires  spaced 
far  apart  offered  large  inductance  (free  space  side).  Using  the  assumption  of  local 
periodicity,  the  inductive  properties  at  any  given  point  were  determined  on  the  tapered 
parallel  wire  grid  using  a  infmite  parallel  wire  grid  model.  On  the  free  space  edge,  there 
was  a  limit  to  the  spacing,  due  to  grating  lobes.  Therefore,  for  near  grazing  incidence  the 
spacing  of  wires  was  kept  to  less  than  X/2  for  the  highest  frequency  of  interest. 

Wait's  Solution  for  Parallel  Wire  Grids.  The  first  analytical  solution  was 
developed  by  Wait  [5].  This  solution  is  for  arbitrary  incidence  on  a  parallel  wire  grid.  In 
his  paper.  Wait  described  the  scattered  electric  field  as  observed  from  a  point  (x,y,z)  as 


(2.1) 


For  a  description  of  a,  d,  0  and  (|>,  sec  Figure  3.  If  the  observation  point  is  in  the  far  zone 
and  assuming  no  grating  lobes,  only  the  m  =  0  term  will  be  of  consequence  (the  higher 
order  terms  will  act  as  evanescent  waves).  This  reduces  the  above  equation  into 


E,,n.iocose^ 

2dcos(|) 

where 

q  =  -  tx|cos(j)cos0  +  ysin())cos0  -  zsin0. 


In  this  particular  study  0=0,  and  therefore  E|and  Ey  will  be  zero.  This  is  reasonable  as  no 
cross  polarization  terms  were  expected  when  scanning  wires  about  the  orthogonal  axis.  In 
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the  above  equation,  q  represents  the  phase  terms  associated  with  propagation  from  the  grid 
to  the  observation  point.  The  current  term  lo  can  be  represented  by 


This  solution  is  valid  for  wires  spaced  far  apart  (d  >  5a).  This  solution  has  been  shown  to 
match  results  from  Henderson's  PMM  code. 


Yoder*s  Solution  for  Strip  Gratings.  Yoder  [6]  developed  another  solution 

to  parallel  strip  gratings  in  his  master’s  thesis.  This  solution  was  based  upon  the  Periodic 
« 

Moment  Method  and  showed  excellent  results  compared  to  PMM  calculations  and  Wait's 
solution.  The  coordinate  axis  used  by  Yoder  was  the  same  system  as  used  in  Henderson’s 
PMM  code  (see  Figure  4).  This  is  the  same  coordinate  system  that  is  used  throughout  the 
rest  of  this  thesis.  The  incident  field  will  be  described  by  two  angles  a  and  q,  where  a  is 
the  angle  that  the  plane  of  incidence  makes  upon  the  xz  plane  (plane  of  the  array)  and  q  is 
the  angle  that  the  incident  E-field  makes  from  the  >y  axis  in  the  plane  of  incidence. 
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Yoder  then  developed  solutions  for  both  parallel  and  orthogonal  scans  of  the  electric 
fields.  The  orthogonal  scans  are  defined  by  a  =:  0®  and  parallel  scans  are  defined  by 
a  =  90®  .  The  two  solutions  developed  by  Yoder  were 


(EJ  = 
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In  this  study,  the  orthogonal  solution  (Equation  2.4)  is  the  only  important  result  for 

the  parallel  wires^  The  parallel  solution  is  important  for  the  thin  strip  grid. 

z 


Figure  4.  PMM  Coordinate  System.  This  is  the  coordinate  system  used  by 
Henderson's  PMM  code  and  Yoder  [2]  [6].  This  coordinate  system  will  be  used 
throughout  this  thesis. 

Primich's  Solution  for  Multiple  Thin  Strip  Gratings.  Primich  [10] 
provided  a  solution  for  multiple  thin  strip  gratings  of  equally  spaced,  identical  strips  in  a 
non-skewed  array.  This  solution  works  for  strip  gratings  that  are  separated  by  free  space 
with  arbitrary  incidence,  Primich’s  solution  works  very  well  and  compares  favorably  to 
experimental  results.  The  first  capacitive  taper  could  be  modeled  using  Primich's  solution, 
but  numerical  techniques  presented  later  will  be  much  easier  to  implement.  The  second 
capacitive  taper  used  multiple  strip  gratings,  but  the  arrays  were  in  skewed  arrangement 
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(the  second  anay  was  offset  from  the  first  array  by  half  a  center-to-center  spacing).  Also 
the  second  capacitive  design  had  a  thin  dielectric  slab  in  between  the  arrays. 

Felsen  and  Carin' s  Solution  for  Truncated  Strip  Gratings.  Felsen  and 
Carin  I7][8][9]  have  authored  several  papers  on  Floquet  mode  adapted  solutions  to 
truncated  periodic  arrays  in  free  space.  Their  solutions  were  based  upon  a  saddle  point 
technique  of  asymptotic  solutions.  These  solutions  were  very  accurate  for  both  near  and 
far  field  solutions.  However,  these  solutions  are  more  complicated  than  Wait  or  Yoder's 
solutions  and  were  not  used  in  this  study. 

Solutions  for  Conductors  on  a  Dielectric  Slab.  There  are  several  solutions 
for  scattering  from  periodic  structures  on  a  dielectric  slab.  Chen  [11]  and  Lee  [12] 

described  scattering  from  periodic  arrays  of  apertures.  These  solutions  use  MoM 

* 

techniques  to  determine  the  scattered  fields.  These  solutions  also  require  that  any  dielectric 
placed  on  the  array  needs  to  be  symmetrically  placed.  Mc.ugomery  [13]  developed  a 
solution  for  arbitrary  placed  dielectric  media.  Montgomery's  solutions  could  model  the 
second  inductive  design,  but  Henderson's  PMM  code  provided  a  much  easier  numerical 
modeling  of  this  design. 

Meanderline  Grids. 

Overview.  Meanderline  grids  were  first  proposed  by  Lemer  [14]  as  circular 
polarizers  due  to  their  special  ability  to  advance  or  retard  and  electric  field  depending  upon 
the  polarization.  Meanderlines  have  the  property  where  the  impedance  is  either  inductive 
or  capacitive  depending  upon  the  polarization  of  the  E  field  (See  Figure  5).  If  the  E  field  is 
parallel  polarized  to  the  grid,  the  grid  appears  as  a  shunt  inductance.  If  the  E  field  is 
orthogonal  polarized  to  the  grid,  the  grid  will  appear  as  a  shunt  capacitance.  This 
difference  causes  a  linearly  polarized  wave,  which  is  polarized  at  45°  to  the  grid,  to  become 
circularly  polarized.  The  incident  linearly  polarized  field  has  equal  components  parallel  and 
orthogonal  to  the  meanderline  that  are  originally  in  phase  with  each  other.  The  component 
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of  the  field  orthogonal  to  the  meanderline  is  advanced  45°  and  the  component  of  the  field 
parallel  to  the  meanderline  is  retarded  45°.  Since  these  components  are  now  equal  in 
magnitude  and  90°  out  of  phase,  a  circularly  polarized  wave  is  emitted.  Several  layers  of 
meanderlines  are  sometimes  used  to  make  the  circular  polarizer  have  better  wideband 
characteristics. 

Meanderlines  Use  in  This  Study.  Meanderlines  have  the  advantage  over 
parallel  wires  to  be  able  increase  the  inductance,  without  having  to  increase  spacing 
between  elements.  Terret  et  al  [  15]  introduced  a  way  of  determining  the  susceptance  using 
MoM  techniques.  Young  et  al  [16]  originally  produced  a  method  of  fmding  analytically  the 
loss  mechanisms  of  meanderline  designs.  Chu  and  Lee  [17]  then  developed  a  full 
analytical  model  of  multi-layered  meanderline  with  arbitrary  incidence.  Skewed  slot  arrays 
were  used  in  this  study  in  lieu  of  meanderlines  for  a  more  continuous  taper  design.  (See 
Figure  6  and  let  S  approach  0) 


Figure  5.  Meanderlines  Impedance  as  a  Function  of  Polarization.  If 
the  E-Field  is  polarized  parallel  to  the  meanderline,  the  meanderline  design  acts  as 
an  inductor.  If  the  E-Field  is  polarized  orthogonal  to  the  meanderline,  the 
meanderline  design  acts  a  capacitor. 
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Figure  6.  Meanderline  Grid.  The  above  figure  defines  the  spacing  of  the 
meanderlfties,  the  periodicity  (pitch)  of  the  design,  the  line  widths  and  the 
alternating  structure  of  the  meanderlines  to  approximate  slots.  If  S  to  approaches  0, 
the  above  design  becomes  a  skewed  slot  array. 

Young's  Solution  for  Insertion  Loss.  Young  [16]  identified  the  loss 
mechanisms  in  meanderline  designs.  He  identified  the  meanderline  loss,  dielectric 
substrate  loss  and  the  loss  of  spacings  between  meanderline  sheets.  The  loss  term 
associated  with  the  meanderlines  is  directly  related  to  the  Q  of  the  inductive  filter.  The  total 


loss  of  the  meanderline  design  is 


ALA  =  -Jf^dB, 

W  Vu.total 


— 1— + 
Qu.total  2\QuJ^ 


where  Qu.totai  is  defined  as 

it  t  \  m 


vu.iotal  uviiiiv\4  oo 

_L_\  +  _J _ + _ 1 _ , 

Qu,c/  Qu  .spacer  Qu, dielectric 


where  (3u,l  is  defined  by 
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ALa  =  total  insertion  loss 
5  s  skin  depth 
p  s  resistivity 
X  =  free  space  wavelength 

In  these  calculations,  Qu,c  is  considered  much  larger  than  Quj^  and  tho-efore  can 
be  ignored.  Also  .the  Q  of  the  spacer  and  dielectric  can  be  ignored  if  the  thickness  is  less 
than  10  mils  (which  is  true  in  this  study).  The  variables  used  above  are  deOned  in  Figure 
6.  If  Dz  was  set  to  H,  the  meanderlines  could  be  modeled  as  a  skewed  slot  array,  which 
will  make  the  design  have  a  more  continuous  taper  and  a  larger  bandwidth.  This  method 
was  not  used  due  to  ease  of  computation  of  Henderson's  PMM  code. 

TerreVs  Solution  for  Scattering  Using  MoM.  Terret  et  al  [15]  developed  a 
moment  method  solution  for  meanderlines.  Hiis  solution  broke  the  individual 
meanderlines  into  6  modes.  Terret's  solution  was  not  chosen  in  this  study  because  the 
PMM  code  proved  to  be  as  accurate  and  computationally  faster. 

Chu  and  Lee*s  Analytical  Solution.  Chu  and  Lee  [17]  offered  an  analytical 
model  for  multi-layered  meanderline  designs.  In  Chu's  solution,  the  effective  sheet 
impedance  of  the  meanderlines  can  be  incorporated  into  a  transmission  line  analogy  for 
both  orthogonal  and  parallel  solutions.  For  the  parallel  polarized  E-field  (which  is  of 
interest  in  this  study)  the  formula  developed  by  Chu  is 


where  K4  and  K5  are  constants  to  be  experimentally  measured. 

Since  the  constants  K4  and  K5  could  not  be  determined,  this  method  was  not  used. 
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Numerical  Solutions  for  Capacitive  and  Inductive  Tapers 

Overview.  This  section  details  the  numerical  solutions  used  in  this  study.  Due  to 
the  complexity  of  the  designs,  only  Wait's  and  Yoder's  methods  may  be  of  use  to 
mathematically  predict  the  inductance  of  parallel  wire  grids  (the  first  inductive  design). 
Therefore,  numerical  methods  offer  the  best  choice  to  model  the  impedance  tapers.  Each 
section  will  outline  the  workings  of  each  method,  the  limitations  of  each  method  and  the 
resources  required  for  each  method.  Each  one  of  these  methods  are  well  accepted  methods 
to  study  electromagnetics.  All  the  numerical  codes  used  have  been  validated  when 
compared  to  measured  data. 

Uniform  Theory  of  Diffraction  (UTD).  Kouyoumjian  and  Pathak  [18] 
developed  the  UTD  in  their  classic  paper  in  1974.  This  solution  used  asymptotic  methods 
to  solve  for  diffraction  from  a  wedge  (see  Figure  7),  without  the  singularities  at  shadow 
boundaries  from  Keller's  Geometrical  Theory  of  Diffraction  (GTD). 


✓ 

Figure  7.  2D  Wedge  Illuminated  by  an  Infinite  Line  Source.  This 
figure  defines  the  values  used  in  the  UTD  calculation.  For  this  study,  WA  =  0  and 

therefore  n  =2.  P  =  5  to  45  X,  and  P  is  very  large  for  far  field  measurements. 
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The  total  field  at  an  observation  point  is  given  by  the  incident  field  plus  the  reflected  and 
diffracted  fields.  The  incident  and  reflected  fields  are  quite  simple  to  find.  Pathak’s 
diffraction  equation  for  both  soft  (orthogonal  to  the  plane  of  incidence)  and  hard  (parallel  to 
the  plane  of  incidence)  polarizations  is 


and  where  the  transition  function  F(x)  is  defined  as 


F(x)  =  2jV5cei»|  e-i^dt,  (2.12) 

it<p)  =  2cos5{?!!^^.  (213) 

where  N*  are  integers  which  satisfy, 

2jcnN*  -  p  =  ±ic,  (2.14) 

where  n  describes  the  wedge  angle  in  radians  by 

n  =  2-S^.  (2.15) 

Because  the  wedge  is  made  of  flat  plates 

L'  =  L™  =  L”*  =  -.  (2. 16) 

p  +  p' 
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Two  codes  for  UTD  solutions  were  written  for  both  hard  and  soft  polarization. 
These  codes  are  included  in  the  Appendix  C.  Note  that  these  two  codes  use  two  subroutine 
(WDC.FOR  and  FTF.FOR)  to  determine  the  wedge  difhracdon  coefficient  and  the  Fresnel 
transition  function,  which  were  developed  at  OSU.  The  exact  values  used  in  the  UTD  code 
are  given  in  Chapter  3. 

Method  of  Moments  (MoM).  MoM  is  one  of  the  most  widely  used  numerical 
methods  in  electromagnetics.  MoM  can  be  used  to  solve  any  linear  operator  equation  by 
treating  the  equation  as  a  system  of  linear  equations.  MoM  determines  the  equivalent 
scattering  currents  and  then  uses  the  electric  or  magnetic  field  integral  equations  to 
determine  the  scattered  fields.  The  currents  can  be  found  by  using  testing  functions  and 
basis  functions,  fn  this  method,  the  equivalent  scattering  currents  representing  small  wires 
or  strips  are  found. 

The  two  codes  were  used  in  this  study  for  the  TMz  (Hz=  0)  and  TEz  (Ez  =  0) 
polarizations.  The  TMz  polarization  code  was  developed  by  myself  and  uses  point  testing 
functions  (delta  functions)  and  pulse  basis  functions  (see  Figure  8).  For  accuracy  the 
length  of  the  modes  was  kept  beneath  O.lX..  This  code  can  handle  PEC  loaded  scatterers  as 
well.  For  this  study,  only  PEC  scatters  were  considered.  The  TEz  code  was  a  singly- 
periodic  code  for  dipoles  [19]. 


Figure  8.  MoM  Basis  Functions.  These  are  graphical  representations  of  the 
two  different  basis  functions  on  a  thin  strip  used  in  this  study. 

This  code  handles  singly  periodic  structures  which  are  periodic  and  infinite  in  the  z 

direction.  This  code  also  can  handle  PEC  and  loaded  structures  as  well  as  active  elements. 
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This  code  used  piece  wise  sinusoid  functions  for  basis  functions,  therefore  for  accuracy  the 
length  of  the  modes  must  be  kept  beneath  O.SX.  (see  Figure  8).  These  codes  required  only 
a  Sun  workstation  to  run  and  each  run  took  around  10  minutes. 

Periodic  Moment  Method  (PMM).  The  periodic  moment  method  is  a 
derivation  of  the  moment  method  which  fmds  the  transmission  and  reflection  coefficients 
for  doubly-infmite  arrays  of  scatterers  using  Roquet's  theorem.  The  code  used  was 
version  3. 1  by  Henderson  [2].  In  this  study  the  code  was  set  up  to  use  piece  wise  sinusoid 
basis  functions.  This  code  can  handle  arrays  of  slots  and  dipoles  as  well  as  slabs  of 
dielectrics.  To  use  this  code,  local  periodicity  was  assumed.  The  code  output  was  the 
plane  wave  reflection  and  transmission  coefficients  for  orthogonal  and  parallel 
polarizations.  This  numerical  method  was  used  to  determine  the  scattering  from  the 
periodic  structures  and  the  results  were  then  used  to  design  the  tapers.  The  code  has  been 
mathematically  and  experimentally  verified  and  was  a  very  easy  code  to  implement 

Finite  Difference  Time  Domain  (FDTD).  The  conceptually  simplest  of  all 
the  numerical  methods  is  FDTD.  FDTD,  first  proposed  by  Yec  [20],  uses  a  central 
differencing  technique  to  solve  for  Maxwell's  equations  directly  in  the  time  domain.  The 
FDTD  solutions  were  quite  memory  and  time  intensive  on  computational  facilities  (the  runs 
required  SO  MB  of  RAM  and  5  hours  of  run  time  on  a  Silicon  Gr^hics  workstation). 
However,  this  method  was  able  to  handle  inhomogenous  media  very  easily  and  the  output 
was  inverse  Fourier  transformed  to  give  frequency  information  from  0  to  70  GHz.  The 
grid  size  was  determined  by  the  minimum  time  step  and  the  smallest  size  of  scatterers  to  be 
modeled.  The  number  of  time  steps  was  determined  by  the  anwunt  of  time  required  for  the 
field  to  completely  propagate  through  the  medium  [21].  The  code  used  was 
TEARHTR.FOR  Version  C  developed  by  Luebbers  [3]  for  the  United  States  Air  Force. 
This  code  used  a  second  order  Liao  boundary  condition  [22]  and  standard  near  to  far  zone 
transform. 
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Summary 

This  chapter  outlined  the  current  knowledge  on  the  subject  of  creating  inductive  ot 
capacitive  periodic  surfaces.  The  analytic  solutions  presented  were  not  used  in  this  study 
due  to  the  computational  efficiency  and  ease  of  numerical  analysis.  The  diffoent  numoical 
analyses  used  in  this  study  were  also  outlined.  The  implementation  of  the  numerical 
methods  will  be  explained  in  more  detail  in  Chapter  3. 
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Chapter  3:  Methodology 


Introduction 

Overview.  This  chapter  describes  the  methodology  used  in  this  study.  In  the 
first  section,  the  taper  design  methodology  is  explained.  The  second  section  covers  how 
the  tapers  were  built  The  third  section  covers  how  the  measurements  were  made.  The  last 
section  covers  the  numerical  analysis  of  the  tapers.  Also  any  assumptions  or  possible 
sources  of  error  are  identified. 

Methodology  in  Design  of  the  Tapers 

Overview.  This  section  describes  how  the  tapers  were  designed.  The  tapers 
designed  had  capacitive  or  inductive  properties.  The  desired  properties  of  a  taper  design 
were  the  impedance  match  both  ends  of  the  taper  (the  ground  plane  edge  and  the  fiee  space 
edge)  and  the  continuous  ^pering  in  the  design  (no  discontinuities).  These  properties  must 
occur  over  an  appreciable  bandwidth  between  2-18  GHz.  The  difficulty  widi  both 
capacitive  and  inductive  designs  was  to  create  either  a  large  enough  ct^acitance  or 
inductance.  Henderson's  PMM  code  [2]  was  used  to  study  the  frequency  range  2-18  GHz 
for  the  crucial  edge  of  the  t^r  (free  space  edge  for  inductive  taper  and  ground  plane  edge 
for  capacitive  taper).  For  each  design,  both  a  long  and  a  short  taper  were  developed.  The 
long  taper  was  2X  at  the  lowest  frequency  (30  cm)  and  the  short  taper  was  IX  at  the  lowest 
frequency  (15  cm).  From  the  Hrst  measurements,  a  second  inductive  and  capacitive  tapor 
were  developed. 

Test  Platform.  Since  the  results  from  this  study  were  to  ultimately  develop  better 
tapering  mechanisms  for  antennas,  the  tapers  needed  to  be  effective  at  near  grazing 
incidence.  To  study  the  effectiveness  of  the  taper  designs,  a  platform  was  tried  to  simulate 
an  semi-infinite  thin  ground  plane  ( a  2D  wed^  where  n  =  2).  The  edge  must  be  at  least 
SX  (75  cm  for  2  GHz)  long  to  simulate  a  2D  edge.  Also  for  a  plane  wave  assumption,  the 
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source  must  be  at  least  SX  away  from  the  edge.  The  AFTT  Model  Shop  built  the  platform 
for  the  measurements.  To  minimize  edge  diffraction  from  the  sides  of  the  ground  plane, 
the  AEL  horn  was  mounted  at  the  apex  of  a  triangular  plate.  The  triangular  plate  was  then 
mounted  on  a  2  inch  styrofoam  substrate  for  stability.  The  rear  of  the  styrofoam  substrate 
had  the  right  angle  edges  cut  ofr  to  form  two  45^  angle  edges  (this  should  reduce  the  small 
amount  of  scattering  from  the  styrofoam  substrate).  The  styrofoam  extended  30  cm 
beyond  the  long  edge  of  the  triangle  to  support  the  impedance  taper  designs.  The  AEL 
horn  antenna  was  then  mounted  by  a  wooden  bracket  to  the  apex  of  the  triangle.  The  entire 
structure  was  then  mounted  a  rotational  platform  by  attaching  a  quarter-inch  plywood  plate 
under  the  styrofoam  substrate.  Since  the  distance  from  the  antenna  to  the  edge  is  more  than 
5X  ,the  plane  wave  incidence  on  the  edge  was  fairly  valid.  Sec  Figure  9  for  the  setup  for 
the  antenna  platform. 

Inductive.  Tapers.  For  the  inductive  tapers,  the  crucial  edge  was  the  free  space 
edge  (high  impedance).  Therefore,  in  the  design  of  inductive  tapo-s,  it  was  important  to 
develop  the  free  space  side  so  that  the  inductance  was  very  high  over  the  frequency  band. 
To  use  Henderson's  PMM  code  to  develop  the  free  space  side  of  the  inductive  taper,  the 
transmission  coefficient  was  computed  for  the  entire  frequency  range  of  interest  (2-18 
GHz).  The  design  for  the  free  space  edge  of  the  taper  should  have  a  high  transmission 
coefficient  (close  to  0  dB)  over  the  bandwidth.  By  determining  the  free  space  side  of  the 
taper,  the  rest  of  the  taper  was  determined  by  linearly  varying  the  taper  variable  (wire 
spacing  in  design  #1  and  slot  length  in  design  #2)  to  the  ground  plane  edge.  The  first 
design  was  made  by  using  the  inductive  proiTerties  of  a  parallel  wire  grid  with  tapered 
spacing.  The  second  design  used  the  inductive  properties  of  thin  slots  by  tapering  slot 
length  in  a  skewed  slot  array. 
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Top  View 


Side  View 


Capacitive  or  Inductive 
Taper 


AEL  Antenna 


A1  Ground  Plane 
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Styrofoam  | - 1 
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Figure  9.  Antenna  Platform  for  Measurements.  This  platform  was 
designed  to  simulate  a  near  grazing  plane  wave  upon  a  2D  edge.  The  Aluminum 
ground  plane  should  be  illuminated  only  at  the  long  edge  (tapered  edge).  The  entire 
structure  was  then  mounted  as  an  antenna  and  measur^  twdi  with  and  without 
tapers  for  both  polarizations.  The  styrofoam  substrate  was  to  provide  stability  for 
the  U^rs  and  ^e  ground  plane.  The  rear  edges  cut  in  the  styrofoam  were  to  reduce 
any  small  amount  of  scattering  from  the  styrofoam  substrate. 
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Inductive  Taper  Design  #/.  A  parallel  wire  grid  with  constant  radius  wires 
with  varying  spacing  was  used  as  the  fust  inductive  taper.  The  largest  spacing  that  had  a 
high  transmission  coefficient  over  the  frequency  domain  was  1  cm  (see  Figure  10).  The 
wires  were  attached  to  a  2  mil  poly-ethylene  plastic  sheet  which  was  then  attached  to  the 
styrofoam  substrate.  The  wires  that  were  available  were  made  out  of  copper  coated  steel 
welding  rods  and  were  91.44  cm  long  and  had  a  1/32  inch  radius.  These  wires  were  then 
attached  to  the  plastic  with  scotch  tape.  The  spacing  was  linearly  increased  from  0  to  1  cm 
over  the  length  of  the  taper  (see  Figure  1 1 ).  A  long  (30  cm)  and  a  short  { 1 5  cm)  parallel 
wire  taper  were  constructed  and  tested. 


Transmission  Coefficent  vs  Frequency  (inductive  Design  #1) 


Frequency  (GHz) 

Figure  10.  Transmission  Coefficient  vs  Frequency  for  Inductive 
Design  #1.  This  is  the  graph  of  the  transmission  coefficient  for  the  free  space 
edge  of  the  first  inductive  taper  computed  using  PMM  at  10°  off  grazing  incidence. 
This  design  is  for  an  parallel  wire  grid  of  constant  radius  (1/32  inch)  wires  with  1 
cm  spacing. 
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Fi^re  11.  Inductive  Taper  Design  #1.  This  was  the  first  inductive  taper 
design.  Each  wire  had  the  same  radius  (1/32  inch).  The  spacing  was  linearly 
increased  from  0  cm  (ground  plane  edge)  to  1  cm  (free  space  end). 

Inductive  Taper  Design  #2.  The  second  inductive  taper  was  developed  by 

using  the  inductive  properties  of  thin  slots.  The  impedance  tap>eT  was  made  by  increasing 

length  of  slots  in  a  skewed  array  over  the  length  of  the  t^r  (See  Figure  13).  The  largest 

slot  that  had  a  high  transmission  coefficient  over  the  frequency  band  was  350  mils  (See 

Figure  12).  This  time  the  design  was  photo  etched  upon  a  10  mil  dielectric  substrate  (8^  = 

4.5).  The  minimum  line  width  for  the  photo  etch  process  was  5  mils;  therefore,  the  length 

of  the  slots  was  linearly  tapered  from  5  mils  to  350  mils.  The  slot  width  was  held  constant 

at  5  ntils.  Two  inductive  tapers  were  made  using  this  design  (15  and  30  cm  long)  and  then 

tested. 
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Transmission  Coefficent  vs  Frequency  (Inductive  Design  #2) 

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19 


Figure  12.  Transmission  Coefficient  vs  Frequency  for  Inductive 
Design  #2.  This  is  the  graph  of  the  transmission  coefficient  for  the  free  space 
edge  of  the  second  inductive  taper  computed  using  PMM  at  10®  off  grazing 
incidence.  This  design  was  a  skewed  slot  array  of  constant  slot  width  (5  mils) 
wires  with  slot  length  350  mils. 

Capacitive  Tapers.  For  the  capacitive  tapers,  the  edge  that  was  difficult  was  the 
ground  plane  edge.  The  difficulty  was  to  develop  a  design  that  had  large  capacitance  (low 
impedance)  over  a  large  bandwidth.  Using  Henderson's  PMM  code,  the  reflection 
coefficient  was  measured  over  the  bandwidth.  The  reflection  coefficient  for  the  ground 
plane  side  should  be  very  close  to  unity  (0  dB)  over  the  bandwidth.  Two  separate 
capacitive  tapers  were  designed.  The  first  capacitive  taper  was  developed  by  having  thin 
strips  of  tapering  lengths  and  constant  center-to-center  spacing.  The  second  capacitive 
taper  was  made  of  two  arrays  of  thin  strips  with  tapered  lengths  and  constant  center-to- 
center  spacing  with  a  thin  dielectric  substrate  slab  in  between  them.  The  second  array  was 
offset  from  the  first  array  by  half  of  a  center-to-center  spacing. 
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Figure  13.  Inductive  Taper  Design  #2.  This  was  the  second  inductive 
taper  design.  Each  slot  had  the  same  width  (5  mils).  The  slot  length  was  linearly 
increased  from  5  mils  (ground  plane  edge)  to  350  mils  (free  space  end). 

Capacitive  Taper  Design  #/.  The  first  capacitive  taper  was  made  by  spraying 

on  thin  strips  using  copper  paint  on  a  thin  sheet  of  paper.  The  strips  had  a  constant  center 

to  center  spacing  with  a  linearly  decreasing  strip  width  on  a  thin  drafting  paper  with  very 

near  unity  dielectric  constant  (See  Figure  14).  The  strips  had  very  high  impedance  due  to 

difficulty  in  painting  small  straight  strips  over  90  cm  long.  Since  these  strips  were  hand 

painted,  the  smallest  gap  between  strips  was  only  0.5  mm.  From  Henderson's  code,  the 

largest  strip  that  had  a  high  reflection  coefficient  was  5.5  mm  (see  Figure  15);  therefore,  the 

center  to  center  spacing  was  6  mm.  The  measured  results  from  this  first  capacitive  taper  are 

very  questionable  due  to  the  crude  manufacturing.  No  short  taper  was  made. 
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Figure  14.  Capacitive  Taper  Deigns  #1  and  #2.  Both  of  the  capacitive 
tapers  had  constant  center-to-center  spacing  between  the  PEC  strips.  The  strip 
widths  were  varied  linearly  from  the  ground  plane  edge  to  the  free  space  edge.  Hie 
first  capacitive  design  strip  width  was  tapered  from  S.S  mm  to  0.5  mm.  The 
second  ctqiacitive  design  strip  width  was  tapored  from  330  mils  to  5  mils.  Also  the 
second  taper  had  another  tapered  array  of  strips  on  die  bottom  of  die  dielectric 
substrate  that  was  offset  by  167.5  mils. 
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Reflection  Coefficent  vs  Frequency  (Capacitive  Design  #1) 

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19 


Figure  15.  Reflection  Coefflcient  vs  Frequency  for  Capacitive 
Design  #1.  This  is  the  graph  of  the  reflation  coefficient  for  the  ground  plane 
edge  of  the  first  capacitive  taper  computed  using  PMM  at  1 0^  off  grazing  incidence. 
This  design  was  a  thin  strip  array  of  constant  center-to-center  spacing  (6  mm)  witfi 
strip  width  5.5  mm. 

Capacitive  Taper  Design  #2.  The  second  taper  was  photo  etched  on  a  10  mil 
dielectric  substrate  (with  a  5  mil  minimum  line  width).  The  dielectric  substrate  had  a 
dielectric  constant  of  4.5  over  2-18  GHz.  The  second  capacitive  taper  design  had  a 
constant  center  to  center  spacing  design  with  varying  strip  widths.  This  design  had  a 
second  array  of  thin  strips  on  the  bottom  of  the  dielectric  substrate  (to  increase  the 
capacitance).  The  bottom  array  was  the  same  as  the  top  array  except  that  it  was  offset  by 
half  of  a  center-to-center  spacing  (See  Figure  14).  The  widest  strip  that  had  a  5  mil  gap  and 
a  large  reflection  coefficient  over  the  entire  bandwidth  was  330  mils  wide  (See  Figure  16). 
The  strip  widths  were  therefore  linearly  tapered  from  330  to  5  mils.  Therefore,  the  center- 
to-center  spacing  for  the  second  capacitive  design  was  335  mils.  For  the  second  capacitive 
design,  two  tapers  were  built  (15  and  30  cm  long)  and  tested  at  OSU-ESL  compact  range. 
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Reflection  Coefficent  vs  Frequency  (Capacitive  Design  #2) 


Frequency  (GHz) 


Figure  16.  Reflection  CoefTicient  vs  Frequency  for  Capacitive 
Design  #2.  This  is  the  graph  of  the  reflection  coefficient  for  the  ground  plane 
edge  of  the  second  capacitive  taper  computed  using  PMM  at  10°  off  grazing 
incidence.  This  design  was  for  two  thin  strip  arrays  of  constant  center-to-center 
spacing  (335  mils),  with  strip  width  tapered  from  330  mils  to  5  mils. 

Experimentally  Measure  the  Data 

Overview.  This  section  explains  how  the  data  was  experimentally  taken.  The 
antenna  platform  and  tapers  were  explained  in  the  earlier  section.  In  each  experiment,  the 
antenna  platform  without  the  edge  taper  was  measured  and  then  the  platform  with  the  t^rs 
were  measured.  The  measurements  were  made  on  two  separate  days.  The  first  day  the 
first  inductive  tapprs  and  the  first  capacitive  taper  were  measured.  On  the  second  day,  the 
second  inductive  and  capacitive  tapers  were  measured.  When  determining  the  effectiveness 
of  the  tapers,  the  comparison  between  with  taper  and  no  taper  measurements  must  be  on  the 
same  day,  because  no  calibration  was  done  between  each  day's  measurements. 
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Compact  Range.  The  compact  range  used  a  parabolic  reflector  to  create  a  far 
field  observation  point  There  are  two  receive  antennas  at  the  focal  point  of  the  reflector, 
one  for  each  principle  polarization.  In  the  center  of  the  range,  there  is  a  platform  with  a 
rotator  which  supports  the  entire  antenna  platform.  (See  Figure  17). 

The  Measured  Data.  The  total  field  was  measured  in  dB  for  a  variety  of  angles 
and  frequencies.  The  angles  recorded  were  the  bistatic  angles,  because  the  source  was 
fixed  for  all  nteasurements  (about  4.7°  from  grazing  off  the  triangular  ground  plane).  The 
first  measurements  were  azimuth  scans  (900-270°  in  0.5°  increments)  at  discrete 
fi’equencies  (2-18  GHz  in  2  GHz  increments).  TTie  second  measurements  were  frequency 
scans  (2-18  GHz  in  20  MHz  increments)  at  discrete  azimuth  angles  (90°,  120°,  and  150°). 
For  the  inductive  t^rs,  the  American  Electronics  Laboratory  (AEL)  source  horn  was 
mounted  to  have  soft  polarization  (E-field  is  parallel  to  the  edge)  on  the  long  edge.  For  the 
capacitive  tapers,  the  AEL  source  horn  was  rotated  to  have  hard  polarization  (E-field  is 
orthogonal  to  the  edge)  on  the  long  edge.  The  hom  was  mounted  on  the  apex  of  the 
ground  plane  by  a  small  wooden  bracket  that  was  epoxied  to  the  structure.  The  antenna 
was  rotated  to  the  proper  polarization  and  then  bolted  into  the  wooden  bracket.  The  taper 
was  attached  to  the  styrofoam  substrate  with  masking  tape.  The  interface  between  the  taper 
and  the  ground  plane  was  taped  over  using  copper  tape  to  reduce  scattering  from  the 
discontinuity.  To  attach  the  first  inductive  design,  the  first  wire  of  the  taper  was  placed 
next  to  the  ground  plane.  To  attach  the  first  capacitive  taper,  the  first  strip  was  taped  to  the 
ground  plane  using  copper  tape.  To  attach  the  second  taper  designs  to  the  ground  plane,  a 
four  inch  thick  copper  strip  was  photo  etched  on  the  dielectric  substrate  before  the  taper 
began.  The  four  inch  thick  copper  strip  was  then  taped  using  copper  tape  to  the  ground 
plane  to  reduce  scattering  from  this  interface. 
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Figure  17.  Compact  Range  Setup  for  Measuring  the  Bistatic 
Diffraction  Pattern.  This  was  the  setup  for  measuring  the  bistatic  diffraction. 
Note  that  the  shadowed  region  was  from  0°  to  180°  and  the  illuminated  region  was 
from  1 80°  to  360°.  The  transmitting  antenna  should  create  an  incident  wave  that  is 
4.7°  off  grazing  incidence  to  the  knife  edge. 


Numerically  Validate  the  Data 

Overview,  Numerical  techniques  were  then  used  to  validate  the  measurements. 
Each  of  the  following  sections  describes  the  numerical  techniques  used  to  validate  the 
measurements.  The  types  of  numerical  techniques  used  were  the  Uniform  Theory  of 
Diffraction  (UTD),  the  Method  of  Moments  (MoM),  the  Periodic  Moment  Method  (PMM), 
and  the  Finite  Difference  Time  Domain  (FDTD)  technique.  For  information  on  the 
limitations  of  each  numerical  technique  see  Chapter  2. 


Uniform  Theory  of  Dtffraction  (UTD).  To  ensure  that  the  results  of  the  no 
taper  nieasureinents  were  valid,  UTD  calculations  for  a  knife  edge  (n  =  2)  were  confuted 
for  comparison.  The  bistatic  data  was  computed  for  both  soft  and  hard  polarizations  for  an 
infinite  line  source  30  inches  (S  to  45  X  for  2-18  GHz)  from  the  knife  edge  at  an  angle  of 
4.7°  from  the  ground  plane.  The  results  were  then  compared  to  the  measured  data  for  the 
no  taper  measurements  to  ensure  that  the  simulation  of  a  knife  edge  is  valid. 

Method  of  Moments  (MoM).  Two  2D  MoM  codes  were  used  to  validate  the 
results  from  the  first  set  of  tapers.  The  first  capacitive  design  was  tested  using  a  code 
developed  by  Skinner  [19].  This  code  used  the  periodic  n[K>ment  method  for  singly  infinite 
arrays  (innnite  in  the  z  direction)  for  incident  fields  that  are  TEz  polarized.  To  simulate  the 
thin  strips,  thin  dipoles  (radius  less  than  X/10)  are  placed  4  radii  apart  from  each  other  in 
the  z  direction,  this  simulated  a  continuous  thin  strip  (See  Figure  18).  Capt  Skinner's 
code  used  piecewise  sinusoid  basis  functions  requiring  each  mode  to  be  no  longer  than  X/2 
.  The  output  of  Skinner's  code  was  the  scattered  field  as  a  function  of  observation  angle. 
The  field  was  incident  at  an  angle  of  170®.  The  observed  fields  were  measured  from  0- 
360®  in  0.5®  increments.  The  antenna  platform  was  modeled  as  a  large  thin  strip  (Im  long) 
with  a  taper  on  one  side.  The  results  from  the  tapered  edge  were  compared  to  the  results 
from  the  large  thin  strip  alone. 

The  first  inductive  taper  was  tested  using  a  code  developed  by  myself  for  the  TMz 
polarization.  The  TMz  code  used  pulse  basis  frinctions  and  point  testing  functions.  To 
model  the  large  thin  strip,  each  mode  was  less  than  AyiO.  To  model  the  thin  wire  taper, 
two  separate  thin  wire  approximations  were  used.  The  wires  that  were  placed  closer 
together  than  5  wire  radii  (5r)  were  represented  by  square  wires  with  side  length  1.45r  [1]. 
For  the  wires  that  were  placed  further  apart  than  5r,  the  wires  were  represented  by  thin 
strips  with  length  4r  [1]  (See  Figure  19).  This  taper  is  described  in  Figure  20.  For  a 
listing  of  the  FORTRAN  code  see  Appendix  D. 


Figure  18.  MoM  Model  for  Capacitive  Design  #1.  This  was  the  MoM 
m^el  used  to  study  the  first  ct^acitive  design.  The  results  from  the  bistatic 
scattering  with  the  taper  were  compared  to  the  results  from  the  bistatic  scattoing 
with  no  taper. 


a  =  1.4Sr  L  =  4r 

Figure  19.  Equivalent  Thin  Wire  Approxinuitions  iced  in  TMz  MoM 
Code.  The  thin  wires  used  in  the  first  inductive  design  were  approximated  as 
square  wires  and  thin  strips  U]. 
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Figure  20.  MoM  Model  for  Inductive  Design  #1.  This  is  the  MoM  model 
us^  to  study  the  first  inductive  design.  The  results  from  the  bistatic  scattering  with 
both  long  and  short  tapers  were  compared  to  the  results  from  the  bistatic  scattering 
with  no  taper. 

Determine  the  Effective  Impedance.  The  next  step  was  to  determine  the 


actual  impedance  taper  by  determining  the  impedance  from  the  reflection  coefficient  This 
solution  translates  the  reflection  coefficient  to  a  thin  sheet  impedance.  Because  the  taper 


designs  are  very  thin  (about  10  mils),  this  approximation  was  valid  up  to  the  first 
resonance,  beyond  that  resonance  another  method  would  have  to  be  used.  The  equation  to 
determine  the  effective  sheet  impedance  Zs  of  the  grid  from  the  reflection  coefficient,  R, 


was 


=  where. 

2R 

Zcjj  =  ZoCOSTl 


Z.,  = 


COST) 


(3.1a) 

(3.1b) 

(3.1c) 


Periodic  Moment  Method  (PMM).  Henderson's  PMM  code  was  then  used  to 
determine  the  impedance  as  a  function  as  length  on  the  taper.  By  assuming  local 
periodicity,  the  impedance  at  any  point  on  the  taper  can  be  approximated  by  using  the 
reflection  coefficients  from  PMM.  For  PMM  code  listings  see  Appendix  A. 

For  the  first  inductive  design,  the  effective  Cr  of  the  array  was  set  near  zero  so  that 
the  shape  of  the  PWS  approximate  overlapping  triangular  functions.  Overlapping 
triangular  functions  are  better  than  the  standard  PWS  for  the  parallel  wires,  because  the 
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current  on  the  wires  should  be  constant  in  the  z  direction  for  plane  wave  incidence.  Dx  was 
set  to  the  spacing  between  wires  (1  cm  at  the  free  space  end).  Dz  was  set  to  0.25  cm  and 
the  wire  radius  was  set  to  1/32  inch.  The  length  to  width  ratio  was  kqpt  greater  than  8: 1  for 
accuracy  of  PMM.  The  incidence  angle  was  defined  by  a  =  0  and  T|  =  80.  (See  Figure  21) 


Dipole  Length 
Dz 


Figure  21.  PMM  Model  for  Inductive  Design  #1.  This  was  the  model  for 
the  parallel  wire  grid  with  variable  spacing.  The  incident  angle  was  defined  by  a  = 
0  and  T]  =  80.  The  overlapping  dipoles  had  PWS  basis  functions.  To  ensure  that 
the  PWS  approximate  triangular  functions,  the  effective  e  was  set  to  near  zero. 


The  PMM  code  was  then  used  to  model  the  first  capacitive  taper.  The  first 
capacitive  taper  (thin  strips)  were  modeled  by  thin  dipoles  placed  4  wire  radii  apart .  The 
wire  radii  of  the  dipoles  were  set  very  small  (0.025  mils)  so  that  the  length  to  width  ratio 
could  be  kept  greater  than  8: 1  (for  PMM  accuracy).  Dx  was  set  equal  to  the  equivalent  wire 
width  of  the  wire  (0. 1  mil),  this  allowed  the  dipoles  to  approximate  a  continuous  thin  strip. 
Dz  was  set  to  the  center-to-center  spacing  (6  mm  at  the  ground  plane  edge).  The  current  on 
each  dipole  was  modeled  using  five  overlapping  PWS  modes.  Since  these  dipoles  were 
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finite  in  extent,  the  Cr  was  not  set  near  zero  because  PWS  offered  the  best  basis  function 

model.  The  incidence  angle  was  defined  by  a  =  90  and  ti  =  80.  (Sec  Figure  22) 

Top  View 

Overlapping  PWS  y  ® 


r/m^m 


Figure  22.  PMM  Model  for  Capacitive  Design  #1.  This  was  the  model 
for  the  thin  strip  grid  with  constant  center-to-center  spacing  and  varying  length. 

The  incident  angle  was  defined  by  a  =  90  and  t)  =  80. 

The  second  inductive  taper  was  then  naodeled  using  PMM.  The  skewed  slot  array 
kept  a  constant  slot  width  (5  mils).  The  Dz  was  set  equal  to  the  length  of  the  slot  (350  mils 
at  the  free  space  end)  plus  the  line  width  (S  mils).  The  Dx  was  set  equal  to  the  slot  width 
plus  the  line  width.  The  skewed  array  was  modeled  upon  a  thin  (10  mils)  dielectric  slab  (e, 
=  4.5).  The  current  on  each  of  the  slots  was  modeled  using  five  overlapping  PWS  modes. 
The  incident  angle  was  defined  by  a  =  90  and  q  =  80.  (Sec  Figure  23) 

The  second  capacitive  taper  was  then  modeled  using  PMM.  The  thin  strips  were 
modeled  as  dipoles  placed  4  wire  radii  (0.025  mils)  ^art  Dx  was  set  equal  to  the 
equivalent  width  of  the  dipoles  (0. 1  mils).  Dz  was  set  equal  to  the  center-to-center  spacing 
(350  mils  at  the  ground  plane  edge)  plus  the  minimum  gap  size  (5  mils).  Next  a  thin 
dielectric  slab  (10  mils)  was  placed  under  the  first  array  (Er  =  4.5).  A  second  identical  thin 
strip  array  was  placed  under  the  dielectric  and  offset  by  half  of  Dz  (167.5  mils).  The 
current  on  the  dipoles  was  modeled  using  five  overlapping  PWS.  The  incident  angle  was 
defined  by  a  =  90  and  q  =  80.  (See  Figure  24) 
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Figure  24.  PMM  Model  for  Capacitive  Design  #2.  This  was  the  model 
for  the  second  capacitive  design.  The  incident  angle  was  defined  by  a  =  90  and  r\ 
=  80. 
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Finite  Difference  Time  Domain  (FDTD).  Lucbber's  FDTD  code  [3]  was 
then  used  to  validate  the  second  capacitive  and  inductive  tapers.  This  code  was  run  on  a 
Silicon  Graphics  machine.  To  ensure  that  the  dielectric  slab  (lO  mils  thick)  could  be 
modeled  by  the  code  the  cell  size  was  set  to  0.4244  mm  by  0.4244  mm  which  gave  time 
steps  of  1  psec.  The  dielectric  constant  for  the  styrofoam  was  set  to  1.02  and  the  dielectric 
constant  for  the  photo  etch  substrate  was  set  to  4.5.  The  entire  grid  size  for  the  FDTD  run 
was  4300  by  700  (the  Liao  boundary  condition  is  generally  considered  valid  for  grids 
which  have  less  than  a  10:1  rectangularity).  The  computer  runs  for  this  model  required 
over  50  MB  of  RAM  and  several  hours  to  run.  The  observation  angles  measured  were 
90°,  120®,  150®,  180®,  210®,  240®,  and  270®.  The  time  domain  results  were  then 
Fourier  transformed  to  echo  width  using  a  companion  program  to  the  Luebber's  code.  The 


results  from  these  runs  are  given  in  Chapter  4.  (See  Figure  25) 


Figure  25.  FDTD  Model  for  Capacitive  Design  #2.  This  was  the  FDTD 
m^el  for  the  second  capacitive  design.  The  cell  size  of  the  grid  was  0.4244  mm 
by  0.4244  mm.  The  time  step  was  1  psec.  The  incident  angle  was  170®  and  the 
observation  angles  were  90®,  120®,  150®,  180®,  210®,  240®,  and  270®.  The  large 
strip  was  1 .5  m  wide. 
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Summary 

This  chapter  covered  the  methodology  used  in  this  study.  The  design  methodology 
was  examined  and  explained.  The  experimental  setup  used  at  OSU-ESL  was  detailed.  The 
numerical  models  used  were  explained.  The  variables  and  angle  defuiitions  for  numerical 
and  experimental  were  explained.  The  results  of  the  numerical  and  experimental  analysis 
are  given  in  Chapter  4. 
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Chapter  4:  Results  and  Analysis 


Introduction 

Overview.  This  chapter  reviews  the  results  from  the  measurements  and  the 
numerical  analysis  done  in  this  thesis.  In  the  flrst  section,  the  results  from  the  PMM  code 
are  presented.  In  the  next  two  sections,  the  results  from  the  measurements  of  the  first  and 
second  taper  designs  are  outlined.  The  measurements  and  UTD  validation  are  then  given  in 
the  next  section.  In  the  last  two  sections,  the  results  from  the  MoM  code  work  and  die 
FDTD  code  are  presented.  The  only  design  not  tested  by  numerical  means  was  the  second 
inductive  taper.  In  each  section,  the  sources  of  error  and  preliminary  analysis  of  results  are 
given. 

PMM  Output 

Overview.  This  section  covers  the  data  generated  by  Henderson's  PMM  code. 
Specifically,  the  data  was  computed  as  a  function  of  position  on  the  taper  and  a  function  of 
frequency.  The  method  of  changing  the  PMM  ouq>ut  from  reflection  coefficients  to 
impedances  through  the  first  resonance  was  given  in  Chapter  3  (Equations  3.1  a-c).  For 
additional  results  from  PMM  and  PMM  source  code  listings,  see  Appendix  A. 

Impedance  Taper  Edges  as  a  Function  of  Frequency.  Each  impedance 
design  had  a  crucial  side  of  the  taper  (free  space  side  for  inductive  taper  and  ground  plane 
side  for  capacitive  taper).  Henderson's  PMM  code  was  used  to  find  the  reflection  and 
transmission  coefficients  as  functions  of  frequencies  for  the  crucial  edges.  The  coefficients 
were  then  transformed  to  impedances.  For  the  inductive  tapers,  the  free  space  edge  needed 
to  have  high  impedance  (high  transmission  coefficient)  over  the  frequency  range.  For 
capacitive  tapers  the  ground  plane  edge  needed  to  have  low  impedance  (high  reflection 
coefficient)  over  the  frequency  range. 
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Modeling  the  termination  edges  made  the  assumption  of  local  periodicity  invalid. 
However,  the  grating  lobe  problems  presented  by  completely  periodic  structures  modeled 
in  PMM,  would  not  necessarily  occur  in  a  tapered  periodic  surface.  Therefore,  the 
assumption  of  local  periodicity  would  predict  a  grating  lobe  condition  in  a  completely 
periodic  structure  that  would  not  occur  in  a  tapered  periodic  structure.  In  predicting  an 
upper  limit  of  design,  the  taper  should  actually  perform  better  than  PMM  predictions. 

The  first  inductive  design's  crucial  edge  was  the  free  space  edge.  The  edge  was 
modeled  using  PMM  as  defined  in  Chapter  3.  The  first  inductive  design  free  space  edge 
had  high  impedance  only  at  15  GHz.  After  IS  GHz  (the  resonance),  the  method  used  to 
determine  the  sheet  impedance  of  the  structure  was  no  longer  valid.  In  the  valid  frequency 
range,  the  real  impedance  was  very  low  in  comparison  to  the  imaginary  impedance,  which 
corresponded  to  a  lossless  impedance.  (See  Hgure  26) 


Impedance  vs  Frequency  (Inductive  Design  «1) 


Fraquency  (GHz) 


Figure  26.  Impedance  vs  Frequency  for  the  Free  Space  Edge  of 
Inductive  Design  #1.  Note  that  only  at  f  =  IS  GHz  does  the  first  inductive 
taper  have  a  high  value  of  impedance. 
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The  second  inductive  design's  crucial  edge  was  the  free  space  edge.  The  edge  was 
modeled  using  PMM  as  defined  in  Chapter  3.  The  edge  displayed  high  impedance  at  8.S 
GHz.  The  bandwidth  of  this  edge  was  much  larger  than  for  the  first  inductive  design, 
meaning  that  the  free  space  edge  of  this  taper  had  high  impedance  over  a  wider  range  of 
frequencies  compared  to  the  first  inductive  taper.  Also  the  real  impedance  was  still  much 
smaller  than  the  imaginary  impedance,  making  the  design  appear  lossless.  The  change  in 
sign  of  the  impedance  at  8.5  GHz  corresponded  to  the  change  from  inductive  to  capacitive 
impedance  at  the  resonant  frequency  of  the  slots.  This  phase  shift  could  cause  time 
dispersion  in  the  time  domain,  by  causing  a  phase  shift  between  frequencies  (in  the  time 
domain  signal)  below  8.5  GHz  and  frequencies  above  8.5  GHz.  The  phase  shift  should 
not  be  a  problem  for  constant  frequency  applications.  (See  Bgure  27) 


Impedance  vs.  Frequency  (Inductive  Design  #2) 
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Figure  27.  Impedance  vs  Frequency  for  the  Free  Space  Edge  of 
Inductive  Design  #2.  Note  that  in  this  design,  the  resonant  frequency  occurred 
around  8.5  GHz. 
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The  first  capacitive  design's  crucial  edge  was  the  ground  plane  edge.  The  edge 
was  modeled  using  PMM  as  defined  in  Chapter  3.  The  edge  had  low  impedance  (less  than 
1(X)  ohms)  over  the  frequency  range  17-20  GHz.  The  taper  sqipeaied  lossless  due  to  very 
low  real  impedance  over  the  entire  frequeiKy  range.  This  taper  was  not  resonant  at  the 
frequencies  of  interest  and  displayed  no  change  in  sign  of  the  impedance.  This  design  was 
constructed  with  a  minimum  line  width  of  O.S  mm.  (See  Figure  28) 


impedance  vs.  Frequency  (Capacitive  Design  #1) 
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Figure  28.  Impedance  vs  Frequency  for  the  Ground  Plane  Edge  of 

Capacitive  Design  #1.  Note  the  high  impedance  at  lower  frequencies. 

The  second  capacitive  design's  crucial  edge  was  the  ground  plane  edge.  The  edge 
was  modeled  using  PMM  as  defined  in  Chapter  3.  The  edge  in  the  second  capacitive 
design  had  much  better  characteristics  than  the  first  capacitive  design.  The  impedance  was 
very  low  Oess  than  5  ohms)  over  the  frequency  range  6-10  GHz.  The  resonance  of  the 
strips  occurred  around  8.5  GHz  and  introduced  a  phase  shift  (which  again  could  be  a 
problem  in  the  time  domain).  The  anti-resonance  at  16  GHz  showed  where  the  ground 
plane  edge  failed,  as  well  as  where  the  impedance  conversion  was  invalid.  (See  Figure  29) 
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Impedance  vs  Fre<|uency  (Capacitive  Design  *2) 


Figure  29.  Impedance  vs  Frequency  for  the  Ground  Plane  Edge  of 
Capacitive  Design  #2.  Note  that  in  this  design,  low  impedance  occulted  in  the 
range  6-10  GHz. 

These  results  showed  how  well  the  crucial  edge  of  each  taper  worked.  One 
important  fact  to  notice  is  that  the  real  part  of  the  impedance  in  each  taper  was  very  low  in 
comparison  to  the  imaginary  impedance  in  all  of  the  tapers.  Therefoie,  the  real  part  of  the 
impedance  will  be  ignored  for  the  remainder  of  this  study.  In  the  first  inductive  t^r,  the 
impedance  of  the  free  space  side  was  low  except  at  15  GHz.  The  impedance  of  the  second 
inductive  design  was  much  higher  than  the  first  inductive  taper  which  should  make  for  a 
better  taper.  The  first  capacitive  taper  had  a  high  impedance  at  the  ground  plane  edge, 
while  the  second  ci^acitive  uq)er  more  closely  matched  the  ground  plane  edge.  The  second 
designs  had  much  better  impedance  characteristics  and  a  much  larger  bandwidth  than  the 
first  designs,  but  the  second  designs  would  cause  dispersion  in  the  time  domain  due  to  the 
phase  shifts. 
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Impedance  and  Reflection  Coefficient  os  o  Function  of  Position  on 
the  Taper.  By  taking  values  along  the  taper  and  assuming  local  periodicity,  the 
impedance  as  a  function  of  position  on  the  taper  was  found  up  to  the  first  resonance.  In 
each  of  these  graphs,  the  reactance  was  graphed  as  a  function  of  length  from  the  ground 
plane  side  to  the  free  space  side.  Due  to  the  8: 1  length  to  width  ratio  required  by  PMM, 
impedances  near  free  space  on  the  capacitive  side  could  not  be  calculated.  Due  to  the 
equivalent  wire  width  used  by  PMM,  impedances  near  the  ground  plane  side  of  the 
inductive  tapers  could  not  be  calculated.  Since  the  conversion  method  for  impedances 
assumed  a  single  resonance,  the  impedances  could  not  be  calculated  for  frequencies  beyond 
than  the  first  resonance.  Also  graphed  are  the  reflection  coefficients  as  a  functicxi  of 
position  on  the  taper,  which  were  valid  over  the  entire  frequency  range.  The  tapers  are 
shoivn  for  frequencies  2, 6, 10, 14,  and  18  GHz.  Additional  frequencies  are  given  in 
Appendix  A. 

In  the  first  inductive  design,  the  spacing  is  varied  from  0  to  1  cm  over  the  length  of 
the  taper.  From  Figure  30,  the  taper  impedance  continuously  increases  from  zero  to  the 
higher  impedance  of  the  taper  edge.  Note  that  for  18  GHz  the  inductive  taper  appeared  to 
end  early,  this  is  due  to  the  grating  lobe  situation  that  occurred  at  frequencies  greater  than 
IS  GHz  preventing  impedance  calculations.  The  high  impedance  side  of  the  taper  over 
most  of  the  frequency  range  was  very  low  in  comparison  to  free  ^ace.  This  low 
impedance  at  the  free  space  edge  indicated  that  this  taper  was  very  poor  fw  frequencies 
below  IS  GHz.  Also  note  that  this  t^r  displayed  no  phase  shift  over  the  entire  taper  in 
the  frequencies  shown.  The  reflection  coefllcient  as  a  function  of  position  on  the  taper  was 
valid  over  the  entire  frequency  range  (see  Rgure  31).  Note  that  the  reflection  coefficient  is 
high  (neai  H  dB)  over  the  entire  Vapcx.  The  higher  frequencies  appeared  to  work  better,  but 
even  the  high  frequencies  did  not  have  smocHh  gradual  tapers  to  free  space,  (see  Hgures 
30  and  31) 
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Figure  30.  Impedance  vs  Length  for  Inductive  Design  #1.  This  has  a 
relatively  low  impedance  over  the  frequency  range. 


Reflection  Coefficient  va  Length  of  Taper  (Inductive  Design  #1) 


Length  of  Taper 

Figure  31.  Reflection  Coeffleient  vs  Length  for  Inductive  Design 
#1.  Note  the  large  discontinuity  at  the  free  space  end. 
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The  first  capacitive  design  had  the  strip  width  varying  from  S.S  mm  to  O.S  mm  over 
the  length  of  the  tap^  with  a  constant  center^to-cenler  spacing  of  6  nun.  From  Figure  32, 
the  taper  did  a  nice  continuous  job  of  tapering  the  impedance  from  zero  to  near  infinity. 
However,  the  ground  plane  edge  had  a  significant  discontinuity  that  is  difficult  to  see  in 
Figure  32.  This  taper  displayed  no  phase  shift  over  the  frequency  range  and  had  high 
impedance  values  for  the  free  space  edge.  The  reflection  coefficient  as  a  function  of 
position  on  the  taper  was  valid  over  the  entire  frequency  range  (see  Figure  33).  The 
discontinuity  between  the  ground  plane  and  the  taper  were  quite  large  (between  12  and  33 
dB).  The  free  space  side  was  much  more  closely  matched  and  should  not  provide  as  much 
scattering  as  from  the  ground  plane  and  taper  interface,  (see  Figures  32  and  33) 


Impedance  ve  Length  of  Taper  (Capacitive  Design  #1) 

Ground  Plane  Free  Space 


Figure  32.  Impedance  vs  Length  for  Capacitive  Design  #1.  This  taper 
had  large  impedances  at  the  free  space  edge. 
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Reflection  Coefficient  vs  Length  of  Taper  (Capacitive  Design  #1) 


Langth  of  Tapor 


Figure  33.  Reflection  Coefficient  vs  Length  for  Capacitive  Design 
#1.  Note  tte  large  discontinuities  at  the  ground  plane  edge. 

The  second  inductive  taper  was  made  by  tapering  slot  lengths  in  a  skewed  array 

from  S  mils  to  350  mils.  Due  to  the  fmite  length  of  the  elements,  resonances  occurred  in 

this  taper.  As  the  frequencies  increased,  the  resonance  slowly  moved  back  (towards  the 

ground  plane)  on  the  taper,  because  the  slot  lengths  decreased  in  that  direction.  At  the 

resonance,  the  impedance  approaches  positive  infinity  and  the  ground  plane  was  exactly 

matched  to  free  space  (see  Figure  34).  On  the  free  space  side  of  the  resonance,  the 

reflection  coefficient  increased  toward  the  free  space  edge  of  the  taper  (see  Figure  35).  The 

taper  worked  best  when  the  resonance  occurred  near  the  end  of  the  taper  (around  8  GHz). 

When  the  resonance  did  not  occur  on  the  taper  (frequencies  less  than  8  GHz),  the  taper  did 

not  work  as  well  because  of  the  high  reflection  coefficient  at  the  free  space  edge.  When  tie 

resonance  occurs  on  the  taper,  but  not  on  the  end,  there  can  additional  scattering  from  the 

taper  free  space  edge  if  the  reflection  coefficient  at  the  free  space  end  was  not  low  enough 

(above  10  GHz).  (See  Figures  34  and  35) 
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Impedance  vs  Length  of  Taper  (Inductive  Design  #2) 
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Figure  34.  Impedance  vs  Length  for  Inductive  Design  #2.  This  taper 
did  have  resonances. 


Reflection  Coefficient  vs  Length  of  Taper  (Inductive  Design  #2) 
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Figure  35.  Reflection  Coefficient  vs  Length  of  Inductive  Design  #2. 
Note  that  only  around  8  GHz  did  the  free  space  si(te  have  a  low  reflection 
coefflcienL 
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The  second  capacitive  design  had  tapered  strip  widths  from  330  mils  to  S  mils.  The 
resonances  occurred  at  approximately  the  same  position  on  the  taper.  The  lower 
frequencies  approached  the  resonance  at  faster  rate  and  have  a  more  gradual  taper  (see 
Fngure  36).  At  the  resonance,  the  impedance  approaches  negative  infinity  and  the  ground 
plane  was  exactly  matched  to  free  space.  On  the  ground  plane  side  of  the  resonance,  the 
reflection  coefficient  decreased  from  zero  to  large  negative  values  (see  Figure  37).  On  the 
free  space  side  of  the  resonance,  the  reflection  coefficient  increased  toward  the  free  space 
edge.  The  reflection  coefficient  was  still  fairiy  low  at  the  end  of  the  taper  over  most  of  the 
frequency  range.  The  ground  plane  interface  was  still  a  scattering  source  for  low 
frequencies  (due  to  low  reflection  coefficients)  and  high  frequencies  (due  to  a  second 
resonance).  This  taper  displayed  the  best  reflection  and  impedance  characteristics  of  all  of 
the  designs  over  the  frequency  range.  (See  Figures  36  and  37) 


Impedance  vs  Length  of  Taper  (Capacitive  Design  #2} 
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Figure  36.  Impedance  vs  Length  for  Capacitive  Design  #2.  This  taper 
does  have  resonances,  but  has  a  higter  bandwidth  than  the  first  capacitive  taper. 
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Reflection  Coefficient  vs  Length  of  Taper  (Capacitive  Design  #2) 

Ground  Plane  Free  Space 


Figure  37.  Reflection  Coefficient  vs  Length  of  Capacitive  Design 
#2.  Note  that  this  taper  has  a  nice  taper  over  most  of  the  frequency  range. 


Measurements  Without  the  Dielectric  Substrate.  Henderson's  PMM  code 
was  then  used  to  determine  the  effect  of  the  dielectric  slab  on  the  second  impedance 
designs.  If  the  dielectric  slab  had  no  effect,  then  some  of  the  mathematical  methods 
mentioned  in  Chapter  2  could  have  been  used  to  check  the  results  from  PMM.  The  results 
from  the  PMM  code  without  the  dielectric  slab  showed  a  marked  difference  in  both  of  the 
second  designs  when  the  dielectric  is  removed  (See  Figures  38  and  39).  Therefore,  the 
mathematical  methods  mentioned  in  Chapter  2  were  not  used  to  solve  the  second  impedance 
designs. 
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Impedance  vs  Frequency  Onductive  Design  #2) 


Figure  38.  Impedance  vs  Frequency  for  Free  Space  Edge  of 
Inductive  Design  #2  Without  the  Dielectric  Substrate.  Note  that  without 
the  dielectric  slab,  the  impedance  changes  dramatically  es"  xially  at  higher 
frequencies. 


Impedance  vs  Frequency  (Capacitive  Design  #2) 


Figure  39.  Impedance  vs  Frequency  for  Ground  Plane  Edge  of 
Capacitive  Design  #2  Without  the  IMelectric  Substrate.  Note  that 
without  the  dielectric  slab,  the  impedance  changes  dramatically. 
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MeasuredResults  on  the  Impedance  Taper  Designs 

Overview.  This  section  covers  the  results  from  the  two  separate  days  of  testing 
at  OSU-ESL.  Since  the  power  ou^ut  of  the  American  Electronics  Laboratory  (AEL)  hom 
was  not  calibrated  each  day,  the  results  from  one  day's  experiment  can  only  be  compared  to 
the  same  day's  data.  Each  pattern  without  a  taper  was  compared  to  UTD  solutions  to  test 
the  ability  of  the  antenna  platform  to  simulate  an  infuiite  2D  knife  edge.  To  smooth  the 
data,  a  moving  median  was  taken  every  7  data  points  (1 .5^  on  each  side  of  the  data  point) 
for  both  the  measured  and  computed  data. 

First  Inductive  Taper.  The  first  inductive  taper  was  measured  and  showed  no 
improvement.  There  was  no  difference  in  the  total  field  pattern  except  at  18  GHz.  The 
improvement  was  slight  at  18  GHz  and  only  in  the  shallow  shadow  region  (less  than  30^  in 
the  shadow  region).  This  improvement  did  correspond  to  the  expectation  from  Figure  31 
that  18  GHz  would  have  the  best  results.  For  the  other  frequencies,  the  wires  used  were 
not  able  to  produce  high  impedance  and  diffraction  reduction  was  expected  to  be  slight 
The  rest  of  the  azimuth  and  frequency  scans  are  in  Appendix  B. 

First  Capacitive  Taper.  The  measurements  made  on  the  first  capacitive  taper 
design  showed  a  very  limited  effect  on  the  scattered  field.  The  tapers  were  crudely  made 
and  therefore  the  results  from  these  tapers  are  very  questionable.  Since  the  impedance  of 
the  strips  was  very  high  (some  were  open  circuits),  the  PEC  assumption  no  longer  held 
true  for  these  tapers.  However,  for  18  GHz  there  was  a  very  slight  reduction  in  the  deep 
shadow  region.  This  reduction  was  less  than  5  dB  and  only  appeared  present  at  18  GHz, 
which  corresponded  to  the  results  ir-esent  in  Figure  33.  The  azimuth  and  frequency  scans 
for  the  first  capacitive  taper  are  in  Appendix  B. 

Second  Inductive  Taper.  The  second  inductive  taper  was  measured.  The 
results  of  this  taper  were  slightly  better  than  the  first  inductive  taper.  First  the  soft  edge 
diffraction  from  the  platform  with  no  taper  was  then  compared  to  UTD  results  for  a  knife 
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edge  (incidence  angle  =  4.7^  off  the  ground  plane).  The  comparison  of  the  UTD  and 
measured  results  showed  additional  scattering  sources  than  the  knife  edge  (Figure  40).  The 
scattering  from  the  knife  edge  was  predominant  more  than  30^  into  the  shadowed  region. 
Figure  41  shows  the  azimuth  cut  of  the  second  inductive  design  at  10  GHz.  Limited 
improvement  was  displayed  in  the  region  less  than  30^  in  the  shadow  region  at  10  GHz. 
The  improvement  did  correspond  to  the  results  from  Rgure  35.  Figures  42, 43,  and  44 
give  the  frequency  scans  of  the  second  inductive  design  at  the  observation  angles  90^, 

120®,  and  150®  respectively.  The  only  reduction  in  diffraction  occurred  in  the  near 
shadow  region  (30®  in  the  shadow  region  or  150®  observation  angle).  This  taper  might 
have  worked  better  in  the  region  less  than  30®  in  the  shadow  region,  but  the  additional 
scattering  sources  in  the  shallow  shadow  region  might  have  overshadowed  the 
effectiveness  of  the  taper.  The  power  levels  for  the  deeper  shadow  region  were  very  low 
and  reduction  in  the  far  shadow  region  was  difficult  to  determine.  The  taper  seemed  to  be 
ineffective  at  most  look  angles  and  frequencies.  The  rest  of  the  azimuth  and  frequency 
scans  for  the  second  inductive  taper  are  in  Appendix  C. 

Total  Field  vs  Azmuth  Angle  (Soft  Polarization) 


Figure  40.  Comparison  of  UTD  and  Measured  Results  From  the 
Knife  Edge.  Note  the  additional  scattering  in  the  region  less  than  30®  in  the 
shadow  region. 


60 


Total  FiaM  (dB)  Total  FMd  (dB) 


Total  Field  vs  Azimuth  Angle  (Inductive  Design  #2) 


Azimuth  Anglo  (dogroos) 

Figure  41.  Azimuth  Cut  of  Inductive  Design  #2  at  10  GHz.  Note  that 
this  design  had  a  limited  improvement  in  the  range  15°  to  30®  in  the  shadow 
region. 

Total  Field  vs  Frequency  (inductive  Design  #2) 
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Figure  42.  Frequency  Sweep  for  Inductive  Design  #2  at  90^  in  the 
Shadow  Region.  This  graph  showed  no  improvement  over  the  entire  frequency 
range  at  9(y)  in  the  shadow  region. 
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Total  Field  vs  Frequency  (Inductive  Design  #2) 
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Figure  43.  Frequency  Sweep  for  Inductive  Design  #2  at  120**  in  the 
Shadow  Region.  This  graph  showed  no  improvement  over  the  entire  frequency 
range  at  1200  in  the  shadow  region. 

Total  Held  vs  Frequency  (inductive  Design  #2) 
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Figure  44.  Frequency  Sweep  for  Inductive  Design  #2  at  ISO®  in  the 
Shadow  Region.  This  graph  showed  no  improvement  over  most  of  the 
frequency  range  at  150®  in  die  shadow  region.  Only  the  frequency  range  3-10  GHz 
showed  moderate  improvement 
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Second  Capacitive  Taper.  The  second  capacitive  taper  showed  a  marked 
improvement  over  a  bare  edge.  First  the  hard  edge  diffraction  from  the  platform  with  no 
taper  was  then  compared  to  UTD  results  for  a  knife  edge  (incidence  angle  =  4.7^  off  the 
ground  plane).  The  comparison  of  the  UTD  and  measured  results  showed  additional 
scattering  sources  than  the  knife  edge  (Figure  45)  in  the  region  less  than  30^  into  the 
shadowed  region.  Figure  46  shows  the  azimuth  cut  of  the  second  inductive  design  at  10 
GHz.  Marked  improvement  was  displayed  in  the  region  mote  than  in  the  shadow 
region  for  both  the  long  and  short  tapers.  Figures  47, 48,  and  49  give  the  frequency  scans 
of  the  second  capacitive  design  at  the  observation  angles  90^^,  120^,  and  ISO^  respectively. 
The  second  capacitive  design  seemed  to  work  well  over  a  wide  frequency  range,  S-IS 
GHz,  which  agreed  with  the  results  from  Figure  37.  The  taper  did  not  seem  to  be  effective 
in  the  region  less  than  30^  in  the  shadow  region  possibly  due  to  the  additiorud  scattering 
sources.  The  test  of  the  azimuth  and  frequency  scans  are  in  Appendix  C. 

Total  Field  vs  Azmuth  Angle  (Hard  Polarliation) 


Figure  45.  Comparison  of  UTD  and  Measured  Results  From  the 
Knife  Edge.  The  results  showed  additional  scattering  when  the  observation  was 
less  than  30^  in  the  shadowed  region. 
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Total  Reid  vs  Azimuth  Angle  (Capacitive  Design  «2) 


Azimuth  Anglo  (dogroos) 

Figure  46.  Azimuth  Cut  of  Capacitive  D^ign  #2  at  10  GHz.  No»e  that 
this  design  showed  improvement  mote  than  30^  in  the  shadow  region. 


Total  Reid  vs  Frequency  (Capacitive  Design  *2) 
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Figure  47.  Frequency  Sweep  for  Capacitive  Design  #2  at  90®  in  the 
Shadow  Region.  This  graph  showed  a  10  dB  reduction  over  the  frequency 
range  6-16  GHz  at  90®  in  the  shadow  region. 


64 


Total  Fiald  (dB)  Total  FMd  (dB) 


Total  Raid  va  Frequency  (Capacitive  Design  #2) 
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Figure  48.  Frequency  Sweep  for  Capacitive  Design  #2  at  120^  in  the 
Shadow  Region.  This  graph  showed  a  S  dB  improvement  over  the  frequency 
range  6-16  GHz  at  120^  in  the  shadow  region. 

Total  Field  vs  Frequency  (Capacitive  Design  #2) 
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Figure  49.  Frequency  Sweep  for  Inductive  Design  #2  at  150^  in  the 
Shadow  Region.  This  graph  showed  no  improvement  over  the  entire  frequency 
range  at  150^  in  the  shadow  region. 
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MoM  Results  on  First  Taper  Designs 

First  Inductive  Taper.  The  first  inductive  taper  was  then  modeled  using  a 
MoM  code  developed  by  myself  with  pulse  basis  functions  and  point  testing  functions.  To 
model  the  effectiveness  of  the  taper,  a  long  (Im)  thin  strip  was  modeled  using  the  code 
both  with  one  edge  tapered  and  without  any  edge  tapo's.  The  scattered  fields  from  the  large 
strip  (Im)  were  then  compared  with  modifications  using  short  and  long  tapers.  The  MoM 
measurements  were  bistatic  azimuth  cuts  at  every  0.5^  (incident  field  at  10^  from  grazing). 
The  measurements  were  made  at  2. 10, 18  GHz.  Each  TMz  MoM  computer  run  took  about 
10  minutes  on  a  Sun  workstation.  The  results  at  10  GHz  showed  negligible  improvement 
over  the  bistatic  scans  (See  Figure  50).  The  remainder  of  the  MoM  data  showed  no 
improvement  (See  Appendix  D  for  the  rest  of  the  data). 

Scattered  Field  vs  Azimuth  Angle  (Inductive  Design  *1) 


Azimuth  Angle  (degrees) 

Figure  50.  Comparison  of  MoM  Results  for  Inductive  Design  #1. 
The  graph  above  shows  no  improvement  for  inductive  taper  design  #1  at  10  GHz. 
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First  Capacitiv*  Taper.  The  first  capacitive  taper  was  then  modeled  using  the 
TEz  MoM  code  developed  by  Skinner  [19].  The  comparison  was  made  between  the  long 
thin  strip  (Im)  with  no  edge  taper  and  the  long  thin  strip  with  one  edge  tapered  (30  cm 
taper).  The  MoM  measurements  were  bistatic  azimuth  cuts  at  every  0.5^  (incident  field  at 
10^  off  grazing).  The  measurements  were  made  at  2, 10, 18  GHz.  Each  TEz  MoM  run 
took  about  IS  minutes  on  a  Sun  workstation.  The  calculations  at  10  GHz  showed  an 
improvement  for  the  first  capacitive  taper  (See  Figure  SI).  The  capacitive  taper  seemed  to 
work  especially  well  at  18  GHz.  (See  Appendix  D  for  the  rest  of  the  data) 

The  MoM  results  followed  the  measured  data.  The  first  inductive  taper  did  not 
work  very  well,  because  of  the  inability  to  increase  inductance  enough  with  straight  wires. 
The  TMz  MoM  code  showed  no  positive  effect  on  the  bistatic  scattering  frtxn  the  first 
inductive  taper.  One  possible  source  of  error  in  the  TMz  code  was  the  use  of  the  thin  wire 
approximation  with  the  close  proximity  of  the  wires.  The  first  capacitive  t^r  appeared 
from  the  early  PMM  work  to  be  an  effective  taper  design.  The  TEz  MoM  results  validated 
that  the  first  capacitive  taper  could  woik  as  a  taper  design.  The  MoM  results  could  not  be 
experimentally  validated,  due  to  the  crude  manufacturing  of  the  first  capacitive  taper  making 
the  measurements  dubious. 


67 


Scattered  Field  vs  Azimuth  Angle  (Capacitive  Design  *1) 


Azimuth  Angl*  (dsgrMs) 

Figure  51.  Comparison  of  MoM  Results  for  Capacitive  Design  #1. 

The  graph  above  showed  iinprovemem  for  capacitive  taper  design  #1  at  10  GHz. 

The  taper  seemed  to  be  effective  from  20°  to  140®. 

FDTD  Results  on  Second  Capacitive  Taper  Design 

Second  Capacitive  Taper.  The  second  capacitive  taper  was  modeled  using  the 
FDTD  code  developed  by  Luebbers.  The  model  used  to  study  the  capacitive  taper  is  given 
in  Chapter  3.  The  result  from  the  FDTD  code  for  the  90®  observation  angle  is  given  in 
Figure  48.  The  FDTD  data  did  support  the  measurements  that  showed  a  reduction  in 
diffraction  due  to  the  second  capacitive  taper. 

Also  FDTD  supported  the  measurements  since  no  noticeable  reduction  was  made  in 
diffraction  until  the  observation  was  greater  than  30®  in  the  shadow  region  in  both  the 
FDTD  model  and  the  measurements.  The  FDTD  results  showed  improvement  for  the  long 
uq)er  over  the  short  taper  only  at  lower  frequencies  Oess  than  6  GHz).  This  result  agreed 
with  the  measurements.  The  type  of  taper  (linear,  binomial,  triangular,  etc.)  should  only 
matter  when  the  taper  is  shorter  electrically.  The  FDTD  results  suggested  that  the  second 
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capacitive  taper  was  effective  in  the  range  2-13  GHz.  This  result  tended  to  support  the 
experimental  data,  where  the  second  capacitive  taper  was  effective  in  reducing  deep  shadow 
fields  in  the  frequency  range  6-16  GHz. 

Also  the  FDTD  results  showed  additional  scattering  for  the  no  taper  design  in  the 
shallow  shadow  region  (which  are  noted  in  Figures  16  and  17).  FDTD  was  then  used  to 
try  and  isolate  the  sources  of  additional  scattering.  First  slope  diffraction  was  considered, 
but  since  the  source  was  a  simple  source,  no  slope  diffraction  could  occur.  Second,  the 
scattering  from  the  2  inch  thick  styrofoam  was  considered.  In  the  FDTD  measurements  of 
the  entire  antenna  platform,  the  dielectric  constant  of  the  styrofoam  was  assumed  to  be 
1.02.  However,  the  actual  dielectric  constant  of  the  styrofoam  used  was  too  low  for 
measurement  accuracy.  Therefore,  the  dielectric  constant  could  have  much  higher  than 
1.02  (possibly  as  high  as  1.2).  Therefore,  FDTD  was  used  to  determine  the  scattering 
from  the  2  inch  thick  styrofoam  alone  (without  the  ground  plane  or  taper  design)  as  a 
function  of  the  dielectric  constant  The  values  of  1.02, 1.1  and  1.2  were  used  for  the 
dielectric  constant  The  FDTD  results  from  the  dielectric  substrate  alone  showed 
considerable  scattering  in  the  region  less  than  30^  in  the  shadow  region  (see  Figure  53)  for 
all  three  dielectric  constants.  Figure  S3  shows  the  scattered  field  as  a  function  of 
observation  angle  and  dielectric  constant  and  the  scattered  field  was  above  5  dB  for  most  of 
the  shallow  shadow  region.  The  observation  angles  were  taken  from  45®  to  -45®  in  5® 
increments.  The  cell  size  was  be  made  much  larger  (1.5  mm  by  1.5  mm),  since  there  was 
no  thin  dielectric  substrate  to  model.  The  FDTD  results  also  showed  that  this  additional 
scattering  occurred  over  the  entire  frequency  range  of  interest  (see  Figure  54).  Therefore, 
at  least  part  of  the  additional  scattering  was  due  to  scattering  from  the  styrofoam  substrate. 
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ScattMWd  Fi«ld  (dB/tn)  ScattMod  Fiald  (dB/lm) 


FDTD  Scattered  Field  vs  Frequency 


Figure  52.  Comparison  of  FDTD  Results  for  Capacitive  Design  #2. 
The  gr<^h  above  shows  improved  for  the  sectmd  capacitive  taper  design  up  to  13 
GHz.  llie  incident  field  was  at  10^  and  the  far  field  observation  angle  was  90*^. 

Scattered  Field  vs  Azimuth  Angle  (Styrofoam  Alone) 


Figure  S3.  Scattered  Field  from  the  Styrofoam  Substrate  Alone  as  a 
Function  of  Observation  Angle.  The  graph  above  shows  the  scattering  from 
the  2  inch  thick  styrofoam  alone  for  different  dielectric  constants  at  10  GHz. 
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Scattered  Field  vs  Frequency  (Styrofoam  Alone) 


FrtqtMncy  (GHz) 


Figure  54.  Scattered  Field  from  the  Styrofoam  Substrate  Alone  as  a 
Function  of  Frequency.  The  graph  above  shows  the  scattering  from  tl^  2  inch 
thick  styrofoam  alone  for  different  dielectric  constants  at  a  look  angle  of  90^. 


Summary 

This  chapter  covered  the  results  and  analysis  of  this  study.  The  tapers  were 
designed  using  Henderson's  PMM  code  and  were  tested  at  a  compact  range.  The  tapers 
were  also  modeled  using  UTD,  MoM  and  FDTD  codes.  The  inductive  tapers  did  not  work 
very  well.  The  first  capacitive  taper  was  crudely  made  and  showed  only  limited 
improvement.  The  second  capacitive  taper  worked  very  well  over  a  wide  frequency  range 
in  both  the  measurements  and  numerical  data.  The  FDTD  code  was  also  used  to  try  to 
identify  sources  of  additional  scattering. 
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Chapter  5:  Conclusions  and  Future  Studies 


Introduction 

Overview.  This  chapter  covers  the  conclusions  from  the  results  of  the  study  and 
suggests  future  work  for  periodic  tapers.  The  success  of  all  four  taper  designs  used  in  this 
study  are  covered  in  the  first  section.  The  future  work  section  covets  suggestions  for 
different  taper  designs  and  better  methodology  for  studying  tapered  periodic  surfaces. 

Success  of  the  Taper  Designs 

Inductive  Tapers.  'H”*  inductive  taper  designs  were  not  very  successful.  The 
first  inductive  taper's  experimental  results  showed  very  limited  diffraction  reduction  in  the 
shallow  shadow  region  at  18  GHz.  The  TMz  MoM  analysis  of  the  first  inductive  taper 
showed  the  design  had  a  very  small  diffraction  reduction.  The  large  diameter  wires  in  the 
first  taper  were  a  problem,  since  the  inductance  of  the  free  space  edge  could  not  be 
increased  enough  by  spacing  alone  (due  to  grating  lobe  problems).  The  second  inductive 
taper  showed  more  promise  to  work  as  a  taper  design  from  the  PMM  results.  The 
experimental  measurements  showed  a  larger  reduction  in  the  shallow  shadow  region  for  the 
second  inductive  taper  as  compared  to  the  first  inductive  taper.  However,  the  reduction 
was  still  .slight  and  only  over  a  small  frequency  range  (8-10  GHz).  The  slot  array  was  also 
unable  to  impedance  match  the  free  space  end. 

Capacitive  Tapers.  The  capacitive  tapers  were  more  successful  in  reducing  edge 
diffraction.  The  first  capacitive  taper  measuierncnts  displayed  a  small  leduction  in  the 
scattering,  but  the  crude  manufacturing  made  these  results  highly  suspect  The  TEz  MoM 
code  model  of  this  taper  did  show  a  slight  improvement  for  the  first  capacitive  design.  One 
good  characteristic  of  this  first  capacitive  design  was  that  the  impedance  was  entirely 
capacitive  over  the  entire  frequency  range.  This  characteristic  is  good  for  time  domain 
applications,  due  to  the  fact  that  the  phase  shift  by  the  taper  is  constant  and  will  not  produce 
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time  dispersion.  The  second  capacitive  taper  was  highly  effective  over  a  large  frequency 
bandwidth  (6-16  GHz).  The  long  version  of  the  second  capacitive  taper  did  work  better 
than  the  short  version  of  the  second  capacitive  taper,  but  only  for  low  frequencies.  The 
second  capacitive  taper  had  the  best  match  of  all  the  tapers  for  both  the  ground  plane  and 
free  space  edge:.  The  second  taper  did  have  resonances  which  would  produce  dispersion 
in  a  lime  domain  signal  This  property,  while  unwanted  for  time  domain  applications,  is 
not  terrible  for  continuous  frequency  applications. 

Future  Work 

Overview.  This  section  covers  future  work  suggested  by  the  results  to  further  the 
development  and  understanding  of  periodic  taper  drisign.  There  are  limitless  types  of 
periodic  structures  which  could  be  used  for  their  inductive  or  capacitive  properties. 
Improvements  could  also  have  been  made  in  the  measurements  and  modeling  used  in  this 
thesis.  There  are  also  several  different  applications  which  were  not  studied  in  this  paper. 

Inductive  Tapers.  The  inductive  taper  designs  were  not  successful,  but  future 
designs  might  be  effective.  One  improvement  to  the  first  inductive  design  would  be  to  load 
the  parallel  wire  grid  at  the  free  space  edge.  If  these  wires  were  loaded,  the  free  space  edge 
of  the  taper  would  have  a  higher  impedance  and  match  free  space  better.  A  loaded  taper 
may  still  have  high  power  applications.  The  power  of  the  antenna  should  be  reduced  at  the 
free  space  end  of  the  taper  so  that  a  high  resistive  (or  inductive)  wire  would  not  heat  up  a 
great  amount.  Also,  the  heat  dissipated  in  a  transient  antenna  would  be  reduced,  because 
the  average  power  of  a  transient  antenna  is  much  lower  than  for  a  continuous  wave 
antenna.  The  parallel  wire  grid  would  not  have  a  phase  shift  over  the  frequency  range  and 
would  work  well  for  lime  domain  applications.  Another  improvement  to  the  first  design, 
would  be  to  use  thinner  wires  (which  have  more  inductance  and  can  be  spaced  closer 
together).  To  improve  the  second  inductive  design,  a  lower  dielectric  con.siant  dielectric 
might  improve  the  inductive  capabilities  of  the  taper.  The  thin  dielectric  slab  acted  as  a 
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capacitive  sheet  and  by  reducing  the  dielectric  constant,  the  capacitance  of  the  design  would 
be  reduced.  Another  improvement  might  be  to  increase  the  width  of  the  slots,  which  would 
increase  the  transmission  coefficient  at  the  free  space  end  of  the  taper.  Wider  slots  would 
also  result  in  a  decrease  in  the  bandwidth.  Finally,  for  future  inductive  tapers,  other 
designs  need  to  be  tested.  There  are  several  ways  to  create  inductance  with  periodic 
geometries  and  these  geometries  can  easily  be  tapered  to  create  new  inductive  tapers. 

Capacitive  Tapers.  The  capacitive  tapers  were  successful  in  reducing  edge 
diffraction.  One  improvement  to  the  Arst  capacitive  taper  would  be  to  change  from  constant 
center-to-center  spacing  to  constant  gap  size  (the  gap  size  is  the  length  of  the  gap  between 
strips).  This  change  would  make  for  a  more  continuous  taper  design  and  increase  the 
bandwidth  of  the  taper.  Another  improvement  to  the  flrst  taper  would  be  to  load  the  free 
space  edge  of  the  taper.  This  would  reduce  the  scattering  from  the  free  space  end  of  the 
taper.  Due  to  the  success  of  the  second  capacitive  taper,  this  type  of  design  shows  good 
possibilities.  To  improve  the  second  capacitive  taper,  one  might  try  the  changes  mentioned 
above  for  the  first  capacitive  taper.  Another  taper  improvement  would  be  to  try  different 
periodic  geometries  that  create  capacitive  impedance.  Discrete  capacitive  elements  would  be 
better  for  bandwidth  and  continuous  tapers,  however,  these  discrete  elements  would 
resonate  and  prevent  time  domain  applications.  Therefore,  depending  upon  the  application, 
the  taper  could  be  improved  in  a  variety  of  ways. 

Measurements.  The  measurements  in  this  study  were  bistatic  measurements. 
Monostatic  (or  backscatter)  measurements  would  also  be  of  interest  especially  for  RCS 
reduction.  Monostatic  measurements  would  be  much  easier  to  perform  and  would  not 
require  such  a  large  structure  for  measurements.  Also  monostatic  measurements  would  not 
require  all  the  assumptions  that  were  needed  in  a  bistatic  study.  If  a  similar  antenna 
platform  was  to  be  used  in  a  future  study,  the  dielectric  constant  of  the  styrofoam  must  be 
measured  for  proper  modeling.  For  a  dielectric  constant  so  close  to  unity,  a  cavity 
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resonance  method  of  measuring  the  dielectric  constant  would  be  needed  instead  of  the 
waveguide  method.  Other  future  worit  would  be  to  attach  the  tapers  directly  to  antennas 
much  like  edge  cards.  Another  major  area  of  study  for  these  periodic  tapers  would  be  in 
the  reduction  in  edge  scattering  from  a  fmite  size  radome.  For  better  modeling  of  the 
tapers,  singly-periodic  PMM  codes  could  be  used  instead  of  Henderson's  PMM  code. 

This  would  allow  the  tapers  be  designed  without  having  to  assume  local  periodicity.  To 
improve  the  modeling  of  the  platform  using  FDTD,  the  entire  structure  (including  the 
source  antenna)  should  be  modeled  using  a  3D  total  field  FDTD  code.  This  would  get  rid 
of  the  plane  wave  incidence  approximation  and  the  2D  approximation.  By  modeling  the 
entire  structure,  the  additional  sources  of  scattering  would  be  easier  to  determine.  There  are 
many  different  measurements  that  could  be  of  interest  using  these  periodic  tapers. 

Overview 

This  study  was  a  preliminary  review  of  the  use  of  capacitive  and  inductive  periodic 
structures  as  impedance  tapers.  In  this  study,  effective  tapers  were  designed,  built, 
modeled,  and  measured.  MoM  provided  a  valid  model  for  the  first  capacitive  and  inductive 
designs.  FDTD  provided  a  valid  model  for  the  more  complex  second  capacitive  taper.  The 
design  technique  of  using  Henderson's  PMM  code  and  assuming  local  periodicity  seemed 
to  model  the  impedance  tapers  well.  Future  work  would  include  better  measurement 
techniques  and  better  periodic  taper  designs. 
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APPENDIX  A 


RESULTS  FROM  PMM  OUTPUT 


This  appendix  contains  die  all  the  results  from  PMM.  The  first  results  are  the  PMM 
source  code  listings.  Also  in  this  appendix  are  all  the  impedance  and  reflection  coefficient 
values  calculated  from  the  PMM  codes. 


This  i.<?  the  PMM  code  listing  for  the  first  inductive  taper. 


1 1 ILE  "Inductive  Design  #1" 

PLOTFELE  wirelOO.plt 
FREQUENCY  1,20,.  1 
ETA  =80 
ALPHA  =  0 

!  This  is  the  Erst  inductive  design  composed 

!  of  parallel  wires  of  constant  radius  with  a 

!  va^dng  spacing.  In  this  study  the  radius 
!  is  1/32  of  an  inch  and  the  spacing  is  varied 
!  from  .3  to  1  cm. 

!  the  spacing  is  given  in  centimeters 

ANGLE  ALPHA,  ETA 

SPACING  =  1.00 

DX  =  SPACING 

DZ  =  .25 

SET  MODE  SINE 

NONSKEWED  GRID  DX,  DZ  =  DZ 
!  width  in  centimeters 

RADIUS=2.54/32 
WIDTH  =4*RADIUS 
DIPOLE  ARRAY  WIDTH,  .001,  (.001,0.) 
NODE  10,0,0 
NODE20,04DZ 
NODE  30,0,-DZ 
SEGMENT  20,10,20 
SEGMENT  30,10,30 
MODE  10,2030 
END  ARRAY 
XEQ 
EXIT 
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This  is  the  PMM  code  listing  for  the  first  capacitive  taper. 


TITLE  "Capacitive  Design  #1 " 

PLOTFUJE  slot055.plt 
FREQUENCY  1,20,1 
ETA  =  80 
ALPHA=90 

!  This  is  the  first  capacitive  design  composed  of 

!  parallel  thin  struts  which  have  constant  center  to 

!  center  spacing  and  varying  length.  The  length  is  varied 
!  from  .5  cm  to  .55  cm  and  the  spacing  is  .6  cm. 

I 

!  the  Iragth  is  given  in  cm 

LENGTH=.55 

HALFLENGTH=LENGTH/2. 

QUARTERLENGTH=LENGTH/4. 

THIRDLENGTH=LENGTH/3. 

EIGHTHLENGTH=LENGTH/8. 

SIXLENGTH=LENGTH/6. 

WIREWIDTH=0.001*2.54 

DX=WIREWIDTH 

DZ=.6 

HALFDZsDZ/2. 

ANGLE  ALPHA,  ETA 
SET  MODE  SINE 

NONSKEWED  GRID  DX,  DZ  =  DZ 

DIPOLE  ARRAY  WIREWIDTH,  .001 
NODE  LOJIALFLENGTH 
NODE  2,0,THIRDLENGTH 
NODE  3,0,SKLENGTH 
NODE  4,0,0 

NODE  5,0,-SDCLENGTH 
NODE  6,0,-THIRDLENGTH 
NODE  7,0,-HALFLENGTH 
SEGMENT  1,1,2 
SEGMENT  2,2,3 
SEGMENT  3,3,4 
SEGMENT  4,4,5 
SEGMENT  5,5,6 
SEGMENT  6,6,7 
MODE  2,1,2 
MODE  3,2,3 
MODE  43.4 
MODE  5,4,5 
MODE  63,6 
END  ARRAY 

PMM$ACCURACY=500(X) 

XEQ 

EXIT 
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This  is  the  PMM  code  listing  for  the  second  inductive  taper. 


TITLE  "Inductive  Design  #2" 

PLOTFILE  IND350.PLT 
NOSHOW  ALL 
SHOW  COEFHCIENTS 
FREQUENCY  1,20,.  1 
ETA  =  80 
ALPHA  =  90 

!  This  is  the  second  inductive  design  composed  of  single  skewed  slot 
!  array  on  a  thin  (10  mil)  dielectric  substrate.  The  width  of  the  slots  is  kq>t 

!  to  the  minimum  line  width  (S  mil).  The  length  of  the  slots  is  varied  from 

!  5  mil  to  350  mil.  The  spacing  from  each  slot  to  the  next  is  the  average  loigth 

!  of  the  slots  plus  the  minimum  line  width. 

!  the  length  is  given  in  inches 

LENGTHIN=0.350 

LENGTH=LENGTHIN»2.54 

INMIN  =  0.005 

CMMIN=INMIN*2.54 

WIDTH=0.005*2.54 

HALFLENGTH=LENGTH/2. 

QUARTERLENGTH=LENGTH/4. 

THIRDLENGTH=LENGTH/3. 

EIGHTHLENGTH=LENGTH/8. 

SIXLENGTH=LENGTH/6. 

ANGLE  ALPHA.  ETA 
DX  =  Wtt)TH+CMMIN 
DZ  =  LENGTH4CMMIN 
SET  MODE  SINE 
SKEWED  GRID  DX,  DZ  =  DZ 
SLOT  ARRAY  WIDTH,  .001 
NODE  l.O^ALFLENGTH 
NODE  2,0,THIRDLENGTH 
NODE  3.0,SDCLENGTH 
NODE  4,0,0 

NODE  5,0,-SIXLENGTH 
NODE  6,0,-THIRDLENGTH 
NODE  7.0, -HALFLENGTH 
SEGMENT  1,1,2 
SEGMENT  2,2,3 
SEGMENT  33.4 
SEGMENT  4,43 
SEGMENT  5,5,6 
SEGMENT  6,6,7 
MODE  2,U 
MODE  3,2,3 
MODE  43.4 
MODE  5,43 
MODE  63.6 
END  ARRAY 

SLAB  0.010*2.54,(4.5,0.).(1..0.) 

XEQ 

EXIT 
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This  is  the  PMM  code  listing  for  the  second  capacitive  taper. 


TITLE  "Capacitive  Design  #2" 

PLOTFILE  cap330.plt 
NOSHOW  ALL 
SHOW  COEFnCIENTS 
FREQUENCY  1,20,.  1 
ErA  =  80 
ALPHA=90 

!  This  is  the  second  capacitive  design  composed  of  two  arrays  of  thin 

!  strips  with  constant  center  to  center  ^>acing  with  a  thin  (Khnil)  dielectric 

!  in  tetween  the  arrays.  The  strip  len^  is  var^  from  .005  mil  to  .330  mil. 

!  The  spacing  is  held  constant  to  .335  mil. 

! 

!  The  length  is  given  in  mils. 

LENGTHIN=0.330 

LENGTH=LENGTHIN*2.54 

HALFLENGTH=LENGTH/2. 

QUARTERLENGTH=LENGTH/4. 

THIRDLENGTH=LENGTH/3. 

EIGHTHLENGTH=LENGTH/8. 

SDCLENGTH=LENGTH/6. 

WIREWIDTH=0.001*2.54 

DX=WIREWIDTH 

DZ=(0.335*2.54) 

HALFDZ=DZ/2. 

ANGLE  ALPHA,  ETA 
SET  MODE  SINE 

NONSKEWED  GRID  DX,  DZ  =  DZ 

DIPOLE  ARRAY  WIREWIDTH,  .001 
NODE  l.OJIALFLENGTH 
NODE  2.0.THIRDLENGTH 
NODE  3.0,SKLENGTH 
NODE  4,0,0 

NODE  5,0,-SDCLENGTH 
NODE  6,0,-THIRDLENGTH 
NODE  7,0,-HALFLENGTH 
SEGMENT  1,1,2 
SEGMENT  2,23 
SEGMENT  3,3,4 
SEGMENT  4,43 
SEGMENT  5,5,6 
SEGMENT  6,6,7 
MODE2,U 
MODE 333 
MODE  43,4 
MODE  5,43 
MODE  63.6 
END  ARRAY 
1 

SLAB  0.010*2.54,(4.5,0.),(1..0.) 
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DIPOLE  ARRAY  WIREWIDTH.  .001 
NODE  1.0.HALFLENGTH+HALFDZ 
NODE  2,0.THIRDLENGTH+HALFDZ 
NODE  3,0.SDCLENGTH+HALFDZ 
NODE  4,0 JiALFDZ 
NODE  5.0,-SDCLENGIH+HALFDZ 
NODE  6,0,  THIRDLENGTH+HALFDZ 
NODE  7.0.-HALFLENGTH+HALFDZ 
SEGMENT  1,1,2 
SEGMENT  2,2,3 
SEGMENT  3,3,4 
SEGMENT  4,4,5 
SEGMENT  5,5,6 
SEGMENT  6,6,7 
MODE  2,1,2 
MODE  3,2,3 
MODE  4,3,4 
MODE  5,4,5 
MODE  6,5,6 
END  ARRAY 

PMM$ACXTJRACY=50000 

XEQ 

EXIT 


The  following  are  the  gn^hs  of  the  impedance  and  reflection  coefficients  vs 
position  on  the  taper  for  all  four  of  the  designs  tested  in  this  study.  The  different  series 
represent  the  frequency  for  the  calculations. 
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APPENDIX  B 


FIRST  MEASURED  DATA 


The  first  data  in  this  appendix  are  the  magnitude  and  phase  azimuth  cuts  at 
frequencies  2, 10,  and  18  GHz.  The  second  data  in  this  appendix  are  the  freqtrency 
sweeps  and  time  domain  plots  at  azimuth  angles  90^,  120^,  and  1S0<^.  The  OSU-ESL 
designator  system  used  letters  to  designate  which  taper  was  measured.  The  following 
letters  correspond  to  the  OSU-DATA  for  the  first  impedance  designs. 


a3204—- Capacitive  Taper  #1 
b3204— -No  Taper  for  Hard  I^larizaticm 
c3204 — Long  Inductive  Taper  #1 
d3204— -Short  Inductive  Taper  #1 
e3204— -No  Taper  for  Soft  Polarization 
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MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-8O.e-7O.Q-6O.G-5O.G-4O.0-3O.O-2O.O-lO.O  0.0  -180.-135.-90.-45.  0.  45.  90.  ' 


a3204Qh0200-atrq  AZIM.  07/23/93  09:27 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

361  90.00  0.50  2.0000  ATN=  0 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-8O.9-7O.0-6O.G-5O.6-4O.G-3O.G-2O.O- 10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  1 


■-8Q.G-70.G-60. 


b3204di0200-a.frq  AZIM.  07/23/93  09:43 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

361  90.00  0.50  2.0000  ATN=  0 


90.  120.  150.  180.  210.  240.  270. 


Q3204Qh1800-a.frq  AZIM.  07/23/93  0927 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 


361  90.00  0.50  18.0000  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


90.  120.  150.  180.  210.  240.  270. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-80.0-70.6-60.0-50.0-40.6-30.6-20,0-10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  1 


b3204ohlG00-a.frq  AZIM.  07/23/93  09:43 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

361  90.00  0.50  10.0000  ATN=  0 


ANGLE  IN  DEGREES 


90.  120.  150.  180.  210. 

_| _ I _ I _ \ _ L_ 


90. 


120. 


150. 


180. 


~r~ 

210. 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-80.0-70.0-60.0-50.0-40.0-30.0-20.0-10,0  0.0  -180-135.-90.-45.  0.  45.  90.  1 


b3204ah1800-a.frq 

AZIM.  07/23/93 

09:43 

TARGET 

CAP.  SLOT  GRID 

AVE=  1 

REF=OFF 

361  90.00 

0.50  18.0000 

ATN= 

90.  i: 

o_l - 

20.  150. 

J _ l-_  ■■ 

180. 

-■  1 

■57-  '06-  551^81::^  '  '  -■(TD  O'Ol-O  Oc-O'Oi'-O  u7-00?=0'09-0  01-0'08^- 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-8O.G-7O.O-6O.G-5O.0-4O.G-3O.G- 20.0-10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  ' 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-80.6-70.G-60.G-50.G-40.e-30.»-20.0- 10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  ' 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-8O.G-7O,0-6O.G-5O.0-4O.G-3O.G-2O.O- 10.0  0.0  -180.-135  -90.-45.  0.  45.  90,  ' 


c3204av1800-a.frq  AZIM.  07/23/93  10.17 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

361  90.00  0.50  18.0000  ATN=  0 


"1 - \ - 1 - r  ^ - 

90.  120.  150.  180.  210.  2- 


ANGLE  IN  DEGREES 

90.  120.  150.  180.  210.  2- 

_J _ I _ I _ ^ _ 1 _ 


ANGLE  IN  DEGREES 


d3204av0200-a.frq  AZIM,  07/23/93  10:35 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

361  90.00  0.50  2.0000  ATN=  0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-8O.G-7O.G-6O.G-5O.0-4O.0-3O.&-2O.O-1O,OO.O  -180.-135-90.-45.  0.  45.  90.  ! 


e3204av0200-a.frq  AZIM.  07/23/93  10:49 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OPF 

361  90.00  0.50  2.0000  ATN=  0 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-80.0-70.0-60.0-50.0-40.0-30.0-20.0-10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  ' 


e3204avl000-Q.frq  AZIM.  07/23/93  10:49 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

361  90.00  0.50  10.0000  ATN=  0 

90.  120.  150.  180. 


210. 


240. 


270. 


-Bo.&yo 


e3204Qvl800-a.frq  AZIM.  07/23/93  10:49 
TARGET  CAP.  SLOT  GRID  AVE=  1 

REF^OFF 

361  90.00  0.50  18.0000  ATN=  0 


a3204fh0900-a  FR<ScAZ  07/23/93  09.27 

T.ARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  90.00  ATN=  0 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

).0  -57.0  -38.0  -19.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010, 


a3204fh1200-o 

TARGET 


FR&AZ  07/23/93  09:27 
CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

610.35  120.01  ATN=  0 


801  2002.30  610.35  120.01 

-25.-20.  -15.  -10.  -5. 


10.  15.  20.  25 


-25.-20. -15.  -10.  -5.  0.  5.  10.  15.  20.  25 

TIME  IN  NANOSECS 


10.  12. 


16.  18. 


lO—L 
‘  2 


6.  8.  10.  12.  14.  16. 

FREQUENCY  IN  GHZ 


-95.0  -76.0  -57.0  -38.0  -19.0  0.0  ,  -0.010  -0.005  0,000  0.005  O.OlO 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-95.0  -76.0  -57.0  -38.0  -19.0  0.0  rO.010  -0.005  0.000  0.005  0.010 


a3204fh1500-Q  FR&AZ  07/23/93  09:27 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  150.01  ATN=  0 


25.-20. -15.  -10.  -5.  0.  5.  10.  15.  20.  25 


TIME  IN  NANOSECS 


FREQUENCY  IN  GHZ 


-95.0  -76.0  -57.0  -38.0  -19.0  0.0  -0.010  -0  005.  0  000  0.005  OOlO 


b3204fh0900-a  FR&AZ  07/23/93  09:43 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  90.00  ATN=  0 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


FREQUENCY  IN  GHZ 


b3204fh1205-o  FRScM  07/23/93  09:43 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  120.51  ATN=  0 


C5-25.-20. -15.  -10.  -5.  O.  5.  10.  15.  20.  25 

§  I .  I _ L_ _ I _ I _ L__J _ \ _ I _ 

CD 


LU 

00 

O 

Q_ 

OO 

LU 

CH 


LU 


OO 


Q_ 


dH - 1 - 1 - \ - 1 - 1 - 1 - 1 - 1 - 1 - r 

^-25.-20. -15.  -10.  -5.  0.  5.  10.  15.  20.  25 

TIME  IN  NANOSECS 

2.  4.  6.  8.  10.  12.  14.  16.  18. 


FREQUENCY  IN  GHZ 


-95.0. . -76.Q .  -.S7.Q_.-.38.n  -19.0  0.0  .  -0.010  -0.005  ,0.000  0.005  O.OlO. 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

95,0  -76,0  -57,0  -38,0  -19,0  0,0  t0,010  -0,005  0,000  0,005  0,010 


b3204fh15l0-a  FR&AZ  07/23/93  09;43 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  151.01  ATN=  0 


25.-20.-15.-10.  -5.  O.  5.  10.  15.  20.  25 


TIME  IN  NANOSECS 


-95.0  .  t-76.0  .-57.0  -38.0  -19.0  0.0  .  -0.010  -0  005  0  000  0.005  0.010 


c3204fv0900-a  FR<kAZ  07/23/93  10:17 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  90.00  ATN=  0 


c3204fv1205-a  FR&AZ  07/23/93  10:17 

TARGET  CAP^  SLOT  GRID  AVE-  1 

REF=OFF 

801  2002.30  610,35  120.51  ATN=  0 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-95.0  -76.0  -57.0  -38.0  -19.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010, 


c3204fvl510-a 

TARGET 


FR&AZ  07/23/93  10;17 
CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

610.35  151.01  ATN=  0 


801  2002.30  610.35  151.01 

-25.-20.  -15.  -10.  -5. 


20. 


25.-20.-15.-10.  -5.  0.  5.  10.  15.  20.  25 

TIME  IN  NANOSECS 

2.  4.  6.  8.  10.  12.  14.  16.  18. 


6.  8.  10.  12.  14.  16. 

FREQUENCY  IN  GHZ 


-95.0  -  76.0  -57.0  -.38.0  -19.0  0.0  -0.010  -0.005  .  0.000  .  0.005  O.OlO 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

i.O  -57.0  -38.0  -19.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010 


d3204fv0900-a  FR&AZ  07/23/93  10:35 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  90.00  ATN=  0 


TIME  IN  NANOSECS 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


FREQUENCY  IN  GHZ 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-95.0  -76.0  -57.0  -38.0  -t9.0  0.0  rO.010  -0.005  0.000  0.005  0.010 


d3204fvl200-a  FR&AZ  07/23/93  10.35 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  120.01  ATN=  0 


TIME  IN  NANOSECS 


-38.0  -19.0  0.0  . -0.010  .. -0.005  0.000  0.005  aOlO 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-95.0  -76.0  -57.0  -38.0  -19.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010 


d3204fv1500-a  FR&A2  07/23/93  10:35 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  150.01  ATN=  0 


TIME  IN  NANOSECS 


FREQUENCY  IN  GHZ 


e3204fv0900-a  FR4cAZ  07/23/93  10:49 

TARGET  CAP.  SLOT  GRID  AVE=  1 

REF=OFF 

801  2002.30  610.35  90.00  ATN=  0 


FREQUENCY  IN  GHZ 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

i.O  -57.0  -38.0  -19.0  0.0  rO.010  -0.005  0.000  0.005  0.010 


FREQUENCY  IN  GHZ 


-95.0  -76,0  -57.0  -38.0  -19.0  0.0  .  '-0.01D  -0  005  0.000  0  005  O.OlO 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-95.0  -76.0  -57.0  -38.0  -19.0  0.0  tO.OIO  -0.005  0.000  0.005  0.010 


e3204fv1500-a  FR&AZ  07/23/93  10:49 

target  cap.  slot  grid  AVE=  1 

REF=OFF 

801  2002.30  610.35  150.01  ATN=  0 


TIME  IN  NANOSECS 


AFPENDIX  C 


SECOND  MEASURED  DATA 

This  a|^)endix  contains  the  UTD  codes  and  data  as  well  as  the  second  impedance 
designs  experimental  results.  The  first  set  of  data  is  die  source  codes  for  the  UTD 
calculation  for  both  princpk  polarizations.  The  second  set  of  data  is  the  UTD  calculated 
data  fot  both  principle  polarizations  at  frequencies  2. 10,  and  18  GHz.  The  third  set  of  dtoa 
is  the  azimuth  cuts  of  the  second  impedance  designs  at  frequencies  2, 4, 6, 8, 10, 12, 14, 
16,  and  18  GHz.  The  last  set  of  data  is  the  frequency  sweeps  and  time  domain  response  of 
the  second  impedance  designs  at  azimutii  angtes  90^,  120^^,  and  ISO^. 


The  UTD  source  code  for  sctft  polarization  (electric  line  source)  is  the  fcdlowing: 


PROGRAM  SOFTPOLARIZA‘nON_DIFFRACnON 
C 

C  THIS  PROGRAM  WORKS  FOR  AN  ELECTRIC  LINE  SOURCE  IN  THE 

PRESENCE  OF  A  2D 

C  WEDGE.  N  IS  THE  WEDGE  FACTOR.  RHO  IS  THE  OBSERVATION 
DISTANCE  IN  WAVELENGTHS 

C  RHOP  IS  THE  DISTANCE  FROM  THE  SOURCE  TO  THE  EDGE  IN 

WAVELENGTHS 
C 
C 

COMPLEX  J,UI,UR,UGTD,UUTD,UIO,UI(}R,DS  J>H,U  1  ,U2,U3,Z 
REAL  PHI,PHIP4>HID,PHIPDRHORHOP,DI,S,L,QR,K,NRSBJSB,X 
REALEl,E2E3,ofrset 

c  OPEN(UNIT=21JTLE=*GO_db_sJCLr,STATUS=’UNKNOWN0 

c  OPEN(UNIT=22JTLE=GTD_db_S.XLr,STATUS=’UNKNOWN0 

OPEN(UNIT=24,FILE=*UTD_DB_SJCLl*,STATUSs*UNKNOWbr) 
c  open(unit-25,files'utdnor_s.xir,statuss'unkiiown’) 

J=(0.,1.) 

PI=3. 14592654 

K*2.*PI 

offset=0. 

C 

N=2. 

write(''‘,*)How  far  away  is  the  source  in  wavelengths? ' 


iead(*.*)RHOP 

wiit^*,*)’What  is  the  offset  in  dB? ' 
read(*.*)offset 

RHOP^. 

PHIPD=atan(2.S^)*  18(VPI 
RHOslOOO. 

RSB:slgO.-PHIPD 
ISB=180.-fPHIPD 
PHIP=PHIPD*PI/180 

UIO=EXP(-(J*K*RHOP)ySQRT(RHOP) 

001001=0360 
PHID=270.-y2. 

PHI=PHID*PI/180 

Xl=(RHO*COS(PHI)-RHOP^OS(PHIP))**2 
X=Xl+<RH0*SIN(PHI)-RH0P*SIN(PfflP))**2 
DI=:SQRT(X) 

QR=RHO*RHOP*SIN(PHI+PHIPy(RHO*SIN(PHD+RHOP*SIN(PHIP)) 
L*SQRT((RHOP*SIN(PHIP))**2+(QR-(RHOP^OS(PHIP)))**2) 
S=5QRT((RHO*SIN(PHI))**2+<(RHO^S(PHI))-QR)**2) 

UIQR=EXP(-(J*K*L)ySQRT(L) 

IF  (PHID.GE.ISB)  THEN 
UI=0. 

Gcrroio 

ELSE 

UI=EXP(-(J*K*DDySQRT(DD 
10  ENDIF 

IF  (PHID.GE.RSB)  THEN 
UR=K). 

GOTO  20 
ELSE 

UR=-UIQR*EXP(-(J*K*S))*SQRT(L/(I>S)) 

20  ENDIF 

Z=UIO*EXP(-J*PI/4ySQRT(2*PI*K)*(-.5)*EXP(-J*K*RHOySQRT(RHO) 

UGTD=Z*((l/COS((PHI-PHIPy2)Hl/COS((PHI+PHIPy2))) 

L=RHO*RHOP/(RHO+RHOP) 

CALL  WDC(DS.DH,L,PHID4*HIPD,90.a^ 
UUTD=UIO*DS*EXP(-(J*K*RHO)ySQRT(RHO) 

U1=UI+UR 

U2=UI+UR+UGTD 

U3=UI+UR+UUTD 

El=CABS(Ul)*SQRT(RHOP) 

E2=CABS(U2)*SQRT(RHOP) 


E3=CABS(U3)*SQRT(RHOP) 

if  (el^.O)  that 

eldb»-80 

else 

eldb«20*logl0(el) 

e^if 

if  (e2.eq.O)  then 
e2dbs-80. 

yJiy. 

e2dbs20*logl0(e2) 
end  if 

if  (e3.eq.O)  then 

e3db=-w. 

dse 

e3db=:20*logl0(e3) 
end  if 

c  WRITE(21*)PHID,eidb 
c  WRrrE(22.*)PHID.e2db 

WRITE(24.*)(180.-PHID);/.e3db+offset 
c  write(25,*)PHID.E3 

100  CONTINUE 
STOP 
END 


The  UTD  code  used  for  hard  polarization  was  the  following: 


PROGRAM  HARDPOLARJZAHONJDIFFRACnON 
C 

C  THIS  PROGRAM  WORKS  FOR  AN  MAGNETIC  LINE  SOURCE  IN  THE 

PRESENCE  OF  A  2D 

C  WEDGE  N  IS  THE  WEDGE  FACTOE  RHO  IS  THE  OBSERVATION 
DISTANCE  IN  W  AVELENG 

C  RHOP  IS  THE  DISTANCE  FROM  THE  SOURCE  TO  THE  EDGE  IN 

WAVELENGTHS 

C  this  program  gives  out  UTD  solutions  but  can  also  give  GO  and  GTD  as 

c  wedl 
C 

COMPLEX  J,UI,UR.UGTD,UUTD,UIO.UIQR»DS4)H,U  1  ,U2,U3^ 
REAL  PHI.PHIP,PHID,PHIPDJlHOjmOP,DI,S,L,QR,K.N.RSB.ISB,X 
REAL  Elf  2JE3.offset 

c  OPEN(UNlT=21,FILE=X30  DB_RXLr,STATUS=UNKNOWN0 

c  OPEN(UNIT=22fILE=t3TD_DB_HJCLl*.STATUS=*UNKNOWN0 

OPEN(UNlT=24,FILE=’UTD_DB^H.XLr,STATUS='UNKNOWNT 
c  open(unita2Sfle='utdnor_h.xir,statuS3:'unknown') 

J*(0..1.) 

PI=3.14592654 
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K*2.*PI 

offsetsO. 

C 

N*2. 

write(*,*)'How  far  away  is  the  source  in  wavelengths? ' 
read(*,*)RHOP 

writ^*,*)'What  is  the  offset  in  dB?  * 
read(*,*)ofFsct 
c  RHOP*5. 

PHIPD»atan(2.Sf30.)*  18(VPI 
RHO=1000. 

RSB=180.-PHIPD 

ISB=180>PHIPD 

PH1P=PHIPD*PI/180 

U10=EXP(-(J*K*RHOP))/SQRT(RHOP) 

001001=0360 

PHID=270.-I/2. 

PHI=PHID*PI/180 

Xl=(RHO*COS(PHI)-RHOP*COS(PHIP))**2 
X=X  l+(RHO*SIN(PHI)-RHOP*SIN(PHIP))**2 
DI=SQRT(X) 

QR=RHO*RHOP*SIN(PHI+PHIPy(RHO*SIN(PHI)+RHOP*SlN(PHIP)) 

L=SQRT((RH0P*SIN(PHIP))**2-KQR-(RH0P^0S(PHIP)))**2) 

S=SQRT((RHO*SIN(PHI))**2+((RHOn:OS(PHl))-QR)**2) 

UIQR=EXP(-(J*K*L)ySQRT(L) 

IF  (PHID.GE.ISB)  THEN 
UI=0. 

GOTO  10 
ELSE 

UI=EXP(-(J*K*DI)ySQRT(DI) 

10  END  IF 

IF  (PHID.GE.RSB)  THEN 
UR=0. 

GOTO  20 
ELSE 

UR=UIQR*EXP(-(J*K*S))*SQRT(L/(L+S)) 

20  ENDIF 

Z=UIO*EXP(-J*PI/4ySQRT(2*PI*K)*(-.5)*EXP(-J*K*RHOySQRT(RHO) 

UGTD=Z*((l/COS((PHI.PHIPy2))+(l/COS((PHI+PHn»y2))) 

L=RHO*RHOP/(RHOfRHOP) 

CALL  WDC(DS4)HX.PHID,PHIPD,90.  JO 
UUTD=UIO*DH*EXP(-(J*K»RHO)ySQRT(RHO) 

U1*UI+UR 

U2=UI+UR+UGTD 
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U3=UI+UR+UUTD 


E1=CABS(U1)*SQRT(RH0P) 

E2=CABS(U2)*SQRT(RHOP) 

E3=CABS(U3)*SQRT(RHOP) 

if  (el.eq.O)  then 

eldb=-80 

else 

eldl>»20*logl0(el) 
end  if 

if  (e2.eq.O)  then 

e2db=:-80 

dse 

e2db=20*logl0(e2) 
end  if 

if  (e3.eq.O)  then 
e3db=-80. 

e3db=20*logl0(e3) 
end  if 

c  WRITE(21*)PHID31db 
c  WRrrE(22  ♦)PHID.E2db 

WRITE(24,*)(180.-PHID)/;^3db+offset 
c  write(25.*)PHID.e3 
100  CONTINUE 
STOP 
END 


The  two  subroutines  used  with  these  codes  were  the  following: 


SUBROUTINE  WDC(CDCS,CDCH3  J»HID,PHIPD,BTDJ=N) 

C 

C  UTD  WEDGE  DIFFRACTION  COEFFICIENT 
C 

C  THIS  ROUTINE  IS  IN  BALANIS,  P.848 
C  rr  CALLS  FUNCTION  FTF(X) 

C 

C  INPUTS:  R  (REAL)  =DISTANCE  PARAMETER  LIN  WAVELENGTHS 
C  PHID  (REAL)  =  OBSERV.  ANGLE  PHI  IN  DEGREES 

C  PHIPD  (REAL)  =  INCIDENT  ANGLE  PHI-PRIME  IN  DEGREES 

C  BTD  (REAL)  :=  OBU(^  INdDEKT  ANGLE  BETA  IN  DEGREES 

C  SEE  BALANIS  P.827  FOR  HCTURE  OF  BETA  ANGLE 

C  FN  (REAL)  =  WEDGE  FACTOR  N  (DIMENSIONLESS) 

C  OUTPUTS: 

C  CDCS  (COMPLEX)  =  SOFT  POL.  DIFF.  COEFFICIENT 
C  CDCH  (COMPLEX)  =  HARD  POL.  DIFF.  COEFHCIENT 
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uuuuu 


COMPLEX  D(4).CT.TERMJT(4).FrF.CDCS.CDCH 
DATA  CT.UTPy(-0.05626977.0.05626977).0. 15915494/ 
DATA  Pl,TPLSML/3. 14159265.6.28318531.0.001/ 
DTORsPyi80. 

PHR=PHID*DTOR 

PHPR*PHIPD*DTOR 

SBO=SIN(DTOR*BTD) 

IF(ABS(FN-0.5).LT.SMLXK)TO  10 

DKL*TPI*R 

UFNslTFN 

BETAR=PHR-PHPR 

TERM=CT/(FN*SBO) 

SGN=1. 

1=0 

2 1=1+1 

DN=UFN*(0.5*SGN+UTPI*BETAR) 

N=DN+SIGN(0.5.DN) 

ANG=TPI*FN*N-BETAR 

A=2.*(COS(0.5*ANG)**2) 

X=DKL*A 
Y=PI+SGN*BETAR 
IF(ABS(X).LT.1£-10)  GOTO  3 
FT(I>.FrF(XyTAN(.5*Y*UFN) 

00104 

3  FT(I)=(0..0.) 

IF(ABS(COS(OJ*Y*UFN)).LT.1.E-3XK)TO  4 
FT(D=(1.7725,1.7725)*SIGN(SQRT(DKL).Y) 
FT(D=(FT(D-(0..2,)*DKL*(PI-SGN*/mG))*FN 

4  SGN=-SGN 
D(I)=TERM*FT(I) 

IF(SGN.LT.0.)GOTO2 

BETAR=PHR+PHPR 

IFa.LE.3X30T02 

CDCS=(D(1)+D(2))-(D(3)+D(4)) 

CDCH=(D(1)+D(2))+<D(3)+D(4)) 

RETURN 

10CDCS=(0..0.) 

CDCH=(0..0.) 

RETURN 

END 


COMPLEX  FUNCTION  FTFCXF) 

THIS  IS  THE  FRESNEL  TRANSITION  FUNCTION  USED  IN  UTD. 

THIS  CODE  COMES  FROM  B  ALANIS,  P.  850 
SEE  COMMENTS  THERE 
COMPLEX  FXX(8)  J^(8).C 
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J  DIMENSION  XX(8) 

DATA  XX/.3..5..7,1..1.5.2.3.4..5.5/ 

DATAa/(0..1.y 

DATAFX/(0.5729,0.2677).(0.6768.0.2682),(0.7439,0.2549), 
1(0.8095,0.2322).(0.873,0.  1982).(0.9240.0. 1 577).(0.9658.0. 1073). 
2(0.9797.0.0828)/ 

DATAFXX/(0..0.).(0.5195.0.0025).(0.3355.-0.0665). 

1(0.2187.-0.0757).(0.127.-0.068).(0.0638.-0.0506). 

2(0.0246.-0.0296).(0.0093.-0.0163)/ 

X=ABS(XF) 

IF(X.GT.5.5)GOTO  1 
IF(X.GT.0.3)GOTO  10 

FrF=((  1 .253, 1 .253)*SQRT(X)-(0..2.)*X-0.6667*X*X)*CEXP(CJ*X) 
GOTO20 
10D011N=2.7 

11  IF(X.LTJCX(N))GOTO  12 

12  FrF=FXX(N)*(X-XX(N))+FX(N) 

GOTO  20 

1  FTF=1.-K757(16.*X*X)-0.75)/X/X+CJ*(0.5-157(8.*X*X))/X 
20  IF(XF.GE.O.)RETURN 
FrF=CONJG(FTF) 

RETURN 

END 


The  following  graphs  are  the  results  from  the  previous  UTD  codes  for  frequencies 
2. 10. 18  GHz.  The  first  one  is  for  an  infinite  electric  line  source  (soft  polarization)  and  the 
second  one  is  for  an  infmite  magnetic  line  source  (hard  polarization). 


UTD  Calculations  for  Soft  Polarization 
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UTO  Calculations  for  Hard  Polarization 


128 


The  next  set  of  data  is  the  magnitude  and  phase  azimuth  cuts  at  differmt  angles 
taken  of  the  second  impedance  designs.  The  last  set  of  data  is  the  frequency  sweeps  and 
time  domain  response  of  the  second  impedance  designs.  The  OSU-ESL  designator  system 
used  letters  to  designate  which  taper  was  measured.  The  following  letters  correspond  to 
the  OSU-DATA  taken  on  the  second  impedance  designs. 

The  azimuth  scans  for  the  second  designs  have  the  following  designations: 
a3225— -No  Taper  for  Hard  Polarization 
b3225 — Long  Capacitive  Taper  #2 
c3225 — Short  Capacitive  Taper  #2 
d3225-— No  Taper  for  Soft  Polarization 
e3225-“Long  Inductive  Taper  #2 
f3225— -Short  Inductive  Taper  #2 

The  frequency  scans  for  the  second  designs  have  the  following  designations: 

g3225— -Short  Inductive  Taper  #2 
h3225 — Long  Inductive  Taper  #2 
i3225-— No  Taper  for  Soft  Polwzation 
j3225 — No  Taper  fcM’  Hard  Polarization 
k322S — Long  Capacitive  Taper  #2 
13225— -Short  Capacitive  Taper  #2 
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MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50^0  -40,0  -30.0  -20  0  -10.0  0,0  -180.-135,-00.-45.  0.  h5  00.  ’ 


90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0-50^0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  135.180. 


a31'2Doh0400-a.tar 

TARGET 


AZIM.  08/13/93  10:33 
cap  grid  AVE=  82 

REF=OFF 

361  90.00  .50  4.0023  ATN=  0 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20,0  -10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  135.180. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

60.0  -50,0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135.-50.  -45.  0.  45,  90.  ' 


a3225oh0800-a.tar 

TARGET 


AZIM.  08/13/93  10:37 
cap  grid  AVE=  82 

REF=OFF 

361  90.00  .50  8.0025  ATN=  0 


-60.0  -50.0  -40.0  -30.0  -20  0  -10.0  0.0’  -180.-135.-90.-4  5 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20,0  -10.0  0.0  -180.-135,-90.-45,  0.  45  90,  1 


o3225ahlOOO-a.tQr  AZU  08/13/93  10:39 
TARGET  cap  grid  AVE=  82 

REF=OFF 

361  90.00  .50  10.0023  ATN=  0 


ANGLE  IN  DEGREES 


90.  120.  150.  180.  210.  240.  27C 


ANGLE  IN  DEGREES 


9>-  06-  sa-ostx  xOO  001-  006-  OOi-  oot-  OOS-  009- 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30,0  -20.0  -10.0  0.0  -180-135.-90.-45.  0.  45.  90.  ' 


o3225ohl200-o.tor 

TARGET 


-60.0  -50.0  -40.0  -30.0  -2 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20,0  -10.0  0,0  -180.-135.-90.-45.  0.  45  90.  135.180. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20,0  -10,0  0.0  -180.-135.-90.-45.0,  45.90.135.180. 


a3225ah1600-a.tQr  AZIM.  08/13/93  10;45 
TARGET  cap  grid  AVE=  82 

REF=OFF 

361  90.00  .50  16.0023  ATN=  0 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 


o3225oh1800-o.tar  A2IM.  08/13/93  10:47 
TARGET  cap  grid  AVE=  82 

REF=0FF 

361  90.00  .50  18.0023  ATN=  0 


ANGLE  IN  DEGREES 


-60.0  -50.0  -40,0  -.mo  -20.0  -10  0  0.0'  ^130-135  -90.  -45, 


b3225oh0200-o.tQr  AZIM.  08/13/93  12:01 
TARGET  long  cap  grid  AVE*  82 

REF=OFF 


361  90.00  .50  2.0023  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


90.  120.  150.  180.  210.  240.  270. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50^0  -40,0  -30,0  -20.0  -10.0  0,0  -180,-135.-90,-45  0  45  90.  ' 


b3225ohu400-a.tar  AZIM.  08/13/93  12:03 
TAR(>£T  long  cop  grid  AVE=  82 

REF=OFF 


90.  120.  150.  180.  210.  240. 


ANGLE  IN  DEGREES 


90.  120.  150.  180.  210.  240. 


90.  120.  150.  180.  210.  240. 


ANGLE  IN  DEGREES 


it-- 06-'5u-081:-  -00  0‘0i-  0  03- O'Oe- o  ot- 0  09-  009- 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

a -40,0  -30.0  -20,0  -10,0  0.0  -180-135.-90.-45.  0.  45.  90.  135.180. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135,-90.-45.  0.  45.  90. 


-60.0  -50.0  -40,0  -30.0  -20  0  -10.0  0.0'  -180-1.35, 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20^0  -10.0  0.0  -180.-135,-90.-45.  0.  45.  90.  ' 


b3225ahi:'00-o.tar  AZIM.  08/13/93  12:11 
TARiSET  long  cop  grid  AVE=  82 

REF=OFF 

361  90.00  .50  12.0023  ATN=  0 


90.  120.  150.  180.  210. 


ANGLE  IN  DEGREES 

90.  120.  150.  180.  210. 

_1 _ \ _ I _ I _ ^ _ 


120. 


240 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40^0  -30.0  -20.0  -10.0  0.0  -180-135.-90.-45.  0.  45,  90. 


b3125<jhl400-a.tQr 

TARGET 


AZIM.  08/13/93  12:13 
long  cap  grid  AVE=  82 

REF=OFF 


ANGLE  IN  DEGREES 


-60.0  -50.0  -4 0.0  -30  0  -20.0  -10.0  0.0 •  •  -180.-135.  -90.  -  45. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180-135.-90.-45.  0.  45,  90.  ' 


b32J5ahl600-a.tar  AZIM.  08/13/93  12:15 
TARGET  long  cop  grid  AVE=  82 

REF=OFF 

361  90.00  .50  16.0023  ATN=  0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50  0  -40.0  -30.0  -20.0  -10.0  0.0  -180-135.-90.-45.  0.  45  90.  ' 


b3225oh1800-a.tar 

TARGET 


90.00 


AZIM.  08/13/93  12:17 
long  cop  grid  AVE=  82 

REF=OFF 

.50  18.0023  ATN=  0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20,0  -10.0  0.0  -180.-135,-90.-45,  0.  45,  90. 


c32:'5ojh0200-a.tar  AZIM.  08/13/93  12:25 
TARGET  short  cop  grid  AVE=  82 

REF=OFF 


361  90.00  .50  2.0023  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


c3225oh0400-a.tar  AZIM.  08/13/93  12:27 
TARGET  short  cop  grid  AVE=  82 

REF=OFF 


361  90.00  .50  4.0023  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


c  32  J  !>ah0600-a.  tar 
TARGET 


AZIM.  08/13/93  1 
short  cap  grid  A\ 


361 
90. 


90.00 
120. 


.50  6.0023 
150. 


90.  120.  150. 


180.  210.  240. 

IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180-135-90.-45.  0.  45.  90.  ' 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20^0  -10.0  0.0  -180.-135.-90.-45.  0.  4;.,  90.  135  180. 


c3225oh 1000 -0.  tor 
TARGET 


AZIM.  08/13/93  12:33 
short  cop  grid  AVE=  82 

REr=OFF 

361  90.00  .50  10.0023  ATN=  0 

90.  120.  150.  180.  210. 

J _ \ _ I _ I _ L_ 


90. 

120.  150. 

180.  210. 

ANGLE 

IN  DEGREES 

90. 

120.  150. 

180.  210. 

ANGLE  IN  DEGREES 


MAGNITUDE  II 

-60.0  -50,0  -40.0  -30.0  -2 


-60.0  -50 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50,0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135.-90.-45.  0.  45,  90.  ’ 


c3225oh1400-a.tar 

TARGET 


AZIM.  08/13/93  12:37 
short  cop  grid  AVE=  82 

REF=OFF 


361  90.00  .50  14.0023  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


cJ2:'bohl600-a,tar  AZIM.  08/13/93  12:39 
TARGET  short  cop  grid  AVE=  82 

REF=OFF 

361  90.00  .50  16.0023  ATN=  0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20^0  -10.0  0.0  -180.-135.-90.-45,  0,  45.  90.  ' 


c3225ahl800-o.tor 

TARGET 


AZIM.  08/13/93  12:41 
short  cop  grid  A\/E=  82 

REF=OFF 

361  90.00  .50  18.0023  ATN=  0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50^0  -40,0  -30.0  -20.0  -10.0  0.0  -180.-135.-90,-45.0.  45.  90, 


<J5225aa0200-Q.tar 

TARGET 


AZIM.  08/13/93  13:10 
bare  v  AVE=  82 

REF=OFF 

36!  90.00  .50  2.0023  ATN=  0 


90. 


120. 


150. 


180. 


210. 


24 


270. 


ANGLE  IN  DEGREES 


10. 


ir-  06-  bti-OSi  X  xO'O  O'Ol-  006-  0  02-  O  Of-  0  OS-  0'09z 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.C -30.0  -20.0  -10.0  0.0  -180.-135.-90.-45.  0.  45.  90.  135.180. 


d3  2  L'  I’laoO 4 00 -  Q .  t  or 
TARGET 


AZIM.  08/13/93  13:12 
bare  v  AVE=  82 

REF=OFF 


361  90.00  .50  4.0023  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50^0  -40.0  -30.0  -20,0  -10.0  0.0  -180-135,-90.-45.  0.  45  90.  ' 


d3225oo0600-a.tar 

TARGET 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135.-90.-45.  0.  “15.  90.  ' 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -I80.-135.-90. -45.  0.  45,  90.  : 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50,0  -40^0  -30.0  -20.0  -10,0  0.0  -180-135.-90.-45.  0.  45  90,  ‘ 


'09- 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50  0  -40.0  -30.0  -200  -10,0  0.0  -180.-135.-90,-45.  0.  45.  90.  ’ 


MAGNITUDE  IN  DB  PHASE  IN  DEGPE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180-135.-90.-45.  0.  45.  9' 


d32  2  5aa  1800  -  a.  tar 
TARGET 


AZIM.  08/13/93  13:26 
bare  v  AVE=  82 

REF=OFF 

361  90.00  .50  18.0023  ATN=  0 


ANGLE  IN  DEGREES 


-60  0  -50  0  -40.0  -30.0  -20.0  -10.0  O.O’  • -180-t:'5.-90. -45.  0. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 


e3225aa0400-a.tar  AZIM.  08/13/93  13:37 
TARGET  long  ind.  grid  AVE=  82 

REF=OFF 

361  90.00  .50  4.0023  ATN=  0 


ANGLE  IN  DEGREES 


-60.0  -50.0  -40.0  -30.0  -20.0  -100  0.0'  “180-135.-90.-45. 


MAGNITUDE 

-60.0  -50.0  -40.0  -30.0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30,0  -20^0  -10.0  0.0  -180.-135.-90.-45.  0.  45,  90,  1 


«3225ooU800-a.tar  AZIM.  08/13/93  13;41 


TARGET 


long  ind.  grid  AVE=  82 

REF=OFF 

90.00  .50  8.0023  ATN=  0 


120. 


150. 


180. 


210 


120.  150.  180.  2' 

ANGLE  IN  DEGREES 


120. 


150. 


180. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50,0  -40.0  -30.0  -20^0 -aO  0.0  -180-135.-90,-45.  0.  45  90.  ' 


e  3  2  2  5  oo  1C)00  -  a.  t  or 
TARGET 


AZIM.  08/13/93  13;43 
long  ind.  grid  AVE=  82 

REF=OFF 

361  90.00  .50  10.0023  ATN=  0 


I  ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40,0  -30.0  -20.0  -10.0  0.0  -180,-135.-90.-45,  0.  45.  90.  ' 


-60.0  -50  0  -4.0.0  -30.0  -20,0  -10.0  O.O’  ^180.-135  -90.-45. 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135.-90.-45.  G.  45  90.  ' 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 


e3225aa1GOO-a.tar  AZIM.  08/13/93  13:49 
TARGET  long  ind.  grid  AVE=  82 

REE=OFF 

361  90.00  .50  16.0023  ATN= 


90. 

120. 

150. 

180. 

ANGLE 

IN  D 

90. 

120. 

150. 

180. 

ANGLE  IN  D 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180.-135,-90.-45.  0.  45,  90  i 


e3225ao1800-a.tQr  AZIM.  08/13/93  13:51 
TARGET  long  ind.  grid  AVE=  82 

REF=OFF 


ANGLE  IN  DEGREES 


90.  120.  150.  180.  210.  240.  270 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180-135.-90.-45.  0.  45.  90.  ' 


f3225oa0200-Q.tar  AZIM.  08/13/93  13:59 
TARGET  short  ind  grid  AVE=  82 

REF=OFF 

361  9000  .50  2.0023  ATN=  0 


ANGLE  IN  DEGREES 


-60-0  -.50.0  -40.0  -.30.0  -20.0  -10.0  0.0'  “  180.-135.-90.  -45. 


f3225oQ0400-a.tQr  AZIM.  08/13/93  14:01 
TARGET  short  ind.  grid  AVE=  82 

REF=0FF 

361  90.00  .50  4.0023  ATN=  0 


90.  120.  150.  180.  210.  240.  270. 


90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10,0  0.0  -180.-135.-90.-45.  0.  45  90.  ' 


f3225aaOGOO-Q.tar 

TARGET 


AZIM.  08/13/93  14:03 
short  ind.  grid  AVE=  82 

REF=OFF 

361  90.00  .50  6.0023  ATN=  0 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30,0  -20.0  -10.0  0.0  -180-135.-90.-45,  0.  45.  90,  1 


f3225oa0800-Q.tQr  AZIM.  08/13/93  14:05 
TARGET  short  ind.  grid  AVE=  82 

REF=OFF 


361  90.00  .50  8.0023  ATN=  0 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 

90.  120.  150.  180.  210.  240.  270. 


ANGLE  IN  DEGREES 


•-60.0  -50.0  -40.0  -.30.0 


I  90.  120.  150.  180.  210.  240.  270.' 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20,0  -10,0  0.0  -180,-135-90.-45.  0.  45.  90,  ' 


f3225ool200-a.tQr  AZIM.  08/13/93  14:09 
TARGET  short  ind  grid  AVE=  82 

REF=OFF 

361  90.00  .50  12.0023  ATN=  0 


1 - \ - 1 - 1 - 1 - 1 - 

90.  120.  150.  180.  210.  240.  270 


ANGLE  IN  DEGREES 

90.  120.  150.  180.  210.  240.  270, 


ANGLE  IN  DEGREES 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  -180-135  -90.-45.  0.  45.  90. 


f3225oa1400-a.tar 

TARGET 


AZIM.  08/13/93  14:11 
short  ind.  grid  AVE=  82 

REF=OFF 

361  90.00  .50  14.0023  ATN=  0 


ANGLE  IN  DEGREES 


-60.0  -50.0  -40.0  -30.0 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 

-60.0  -50^0  -40.0  -30.0  -20.0  -aO  0.0  -180-135.-90.-45,  0  45.  90, 


f3225Qal800-a.tQr  AZIM.  08/13/93  14:15 
T ARISE T  short  ind.  grid  A\/E=  82 

REF=OFF 

361  90.00  .50  18.0023  ATN=  0 


ANGLE  IN  DEGREES 


g3225fa0900-a.tor  FREQ.  08/13/93  14:20 
TARGET  freq  s  ind  grid  AVE=  82 

REF=OFF 

1601  2002.30  10.00  90.00  ATN=  0 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


'i 


g3225fa1200-Q.tQr  FREQ.  08/13/93  14:22 

TARGET  freq  s  ind  grid  AVE=  82 

REF=0FF 

1601  2002.30  10.00  120.00  ATN=  0 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-80.G-70.0-60.0-50.&-40.e-30.0-20.0- 10.0  0.0  tO.003-0.002-0.001  0.000  0.001  0,002  0.003 

I  I  I  I  I  I  I  I  I  I  I  I  >  I  I  I  I  ' 


g3225fa1500-a.tar 

TARGET 


FREQ.  08/13/93  14:23 
freq  s  ind  grid  AVE=  82 

REF=OFF 

1601  2002.30  10.00  149.99  ATN=  0 


TIME  IN  NANOSECS 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


1 


h3225fa0900-Q.tar  FREQ.  08/13/93  14:34 
TARGET  freq  I  ind  grid  AVE=  82 

REF=0FF 

1601  2002.30  10.00  90.00  ATN*  0 


o 

o 

o 

o 


o 

oo 

I 


FREQUENCY  IN  GHZ 


h3225fa1500-a.tar  FREQ.  08/13/93  14:36 

TARGET  freq  I  ind  grid  AVE=  82  l 

REF=OFF 

1601  2002.30  10.00  149.99  ATN=  0 


ii225fa0900-a.tar  FREQ.  08/13/93  14:43 
TARGET  freq  naked  v  AVE=  82 

REF=OFF 

1601  2002.30  10.00  90.00  ATN=  0 


kt-20. 

o 

o - 

CD 


-10. 


CNI 

o 

t  .  I  o 
tn  o 
:z:  ^ 

(Tj  Q 

Q-  o 
CT) 

LU  O 

Oc:  8 
o 


s 


Q_  csi 
’S:  O 
O 


-10.  0.  10. 

TIME  IN  NANOSECS 

6.  8.  10.  12.  14. 


20. 


QD  S 
Q  o, 

CNJ 


—  o 
ro 

LU  cb 

q  o 
ZD 

I —  cb 


^  o 


o_ 

cb 

C)__. 

00  1 

I  2 


8.  10.  12.  14. 

FREQUENCY  IN  GHZ 


i3225fal200-a.tor 

target 


FREQ.  08/13/93  14:40 
freq  noked  v  AVE=  82 

REF^OFF 

1601  2002.30  10.00  120.00  ATN=  0 


TIME  IN  NANOSECS 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


i3225fa1500-a.tor  FREQ.  08/13/93  14:41 
target  freq  naked  v  AVE=  82 

REF=OFF 

1601  2002.30  10.00  149.99  ATN=  0 


kt-20.  -10.  0.  10.  20^, 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


t 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010 


j3225fo0900-o.tar 

TARGET 


FREQ.  08/13/93  14:51 
freq  noked  h  AVE=  82 

REF==OFF 

1601  2002.30  10.00  90.00  ATN=  0 


TIME  IN  NANOSECS 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60,0  -50.0  -40.0  -30.0  -20,0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010 


j3225fQl200-a.tar 

TARGET 


FREQ.  08/13/93  14:52 
freq  naked  h  AVE=  82 

REF=OFF 

1601  2002.30  10.00  120.00  ATN=  0 


TIME  IN  NANOSECS 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


FREQUENCY  IN  GHZ 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0,010 


j3225fQl500-a.tor 

TARGET 


FREQ.  08/13/93  14:53 
freq  naked  h  AVE=  82 

REF=OFF 

1601  2002.30  10.00  149.99  ATN=  0 


TIME  IN  NANOSECS 


EREQUENCY  IN  GHZ 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010, 


k3225fa0900-a.tar 

TARGET 


1601  2002.30  II 

-30.  -20. 


FREQ.  08/13/93  14:57 
freq  I  cap  grid  AVE=  82 

REF=OFF 

10.00  90.00  ATN=  0 


-10. 


30. 


20.  -10.  O.  10.  20. 

TIME  IN  NANOSECS 

6.  8.  10.  12.  14.  16, 


2.  4.  6,  8.  10.  12.  14.  16. 

FREQUENCY  IN  GHZ 


-60.G -50.0  -40.Q -30.0  -20.0  -100  00  -0.010  -  -0  005  0.000  0  005  0.010 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010 


k3225fo1200-Q.tim  FREQ.  08/13/93  15:00 
TARGET  freq  I  cap  grid  AVE=  82 

REF=OFF 

1601  2002.30  1220.70  120.00  ATN=  0 


TIME  IN  NANOSECS 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


FREQUENCY  IN  GHZ 


OlU  O  sou  O  '  uOO  O  SOO  O-  OlO'O-  0  0  0  ot- G'Oc- 0  0£- O'O^- 0  OS-  0  09- 


k3225fa1500-a.tar 

TARGET 


FREQ.  08/13/93  14:59 
freq  I  cop  grid  AVE=  82 

REF=OFF 

1601  2002.30  10.00  149.99  ATN=  0 


TIME  IN  NANOSECS 


2.  4.  6.  8.  10.  12.  14.  16.  18. 


i 


-60.0  -50.0  -40.0  -  30.0  -20.0  -10.0  .  0.0  -0.010  -0.005  OOOO  0.005  O.OlO 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0.010, 


I3225fo0900-o.tar 

target 


FREQ.  08/13/93  15:06 
freq  s  cap  grid  AVE*  82 

REF=OFF 

1601  2002.30  10.00  90.00  ATN=  0 


FREQUENCY  IN  GHZ 

i 


I 


I3225fo1200-Q.tar  FREQ.  08/13/93  15:07 
target  freq  s  cop  grid  AVE=  82 

REF=OFF 

1601  2002.30  10.00  120.00  ATN=  0 


MAGNITUDE  IN  DB  IMPULSE  RESPONSE 

-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  rO.OlO  -0.005  0.000  0.005  0,010 


l3225fa1500-a.tQr  FREQ.  08/13/93  15:08 
TARGET  freq  s  cop  grid  AVE*  82 

REF=OFF 

1601  2002.30  10.00  149.99  ATN=  0  I 


30.  -20.  -10.  0.  10.  20.  30 


TIME  IN  NANOSECS 


FREQUENCY  IN  GHZ 


-60.0  -50.0  -40.0  -30.0  :r20.D  -10.0  0.0  .  -0.01C  -aO.05  0.000 .  0.005..  O.OlO 


APPENDIX  D 


METHOD  OP  MOMENTS  DATA 


The  Method  of  Moments  data  was  calculated  by  two  separate  codes.  Thefirstcode 
was  a  TEz  code  written  by  Captain  J.  P.  Skinner  [19].  This  code  was  a  singly-periodk 
PMM  code  which  used  piecewise  sinusoid  basis  functions.  The  second  code  was  a  TMz 
code  written  by  myself  which  used  pulse  basis  functions.  The  fidlowing  is  a  listing  of  the 
TMz  code. 


c 

c  Lt  Russell  A.  Burleson  TEz  2D  scanning  code 

c 

c  Moment  Method  Code  for  2D  geometries  read  from  input  file  TNPUT.DAT 

c  for  TMz  incident  wave, 
c  ida  =  the  number  of  points  aDowed 
c  w  s  firequmcy  in  rad/sec 

c  znraUmat,vinan4mpedanoe,cutient  and  voltaic  matrix 

c  which  ate  sent  out  to  zmatdat  imatdat  vmatdat 

c  zmatdat  can  be  teadin  by  setting  ireadzsl 
c  fteq  =  in  Hz 

c  an  x,y  coordinates  in  meters 

c  o2d.<rat  s  is  the  scattered  echo  width  in  dB  (not  dBm  like  before) 
c  ouQHitdat  =  is  the  scanned  field  in  V/m 
c  inputdat  s  made  from  anothn  program 
c 
c 

parametn  (ida^dlS) 
c 

real  w,uo,pi.phi,phip,phase,degstart,degstop 
teal  x(ida),y(ida).deg,o2d,o2dlog,wavetengtti 
c 

complex  ef(ida),es,sum 

complex  za(ida),zb(ida),zmat(ida,ida),vmat(ida)4mat(ida)J 
c 

integer  nm,itdz,ngo4sc,p,check4wz4wi,iwv 
integer  ia(ida)4b0da) 
c 

open(unit=26,files'o2d.dat',status='unknown') 

open(units27,files'outputdat',status='unknown') 

c 

pi=:3.141592654 

uo^l.256637061e-6 
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j=(0..i.) 

fotmat(3x,lflS.3) 


c 

c 

«* 


C 

«« 

c  THIS  IS  THE  MAIN  PROGRAM 

c 

c 

c  THIS  GETS  THE  INPUT  DATA  FROM  TNPUT.DAT 

C 

C 

call  getinput(wJsc,deg.ngo,irdz,p,phip4mi^,y,ia,ib,za,zb, 

1  ida,iwz,iwi,iwv,degstart,degstop) 
close  (iinit=21) 

C 

writeC*,*)!!!^  is  set  up  to  run  for  '.nm,‘  number  of  modes' 

if  ((nm+O.gtida)  then 

write(6,78) 

78  format(2x,'This  will  not  work-too  many  modes-reset  IDA3 
goto  ^8 
endif 
C 

C  THIS  FELLS  IN  THE  IMPEDANCE  MATRIX  OR  READS  rr  IN  FROM  A  FILE 

C 

C 

if  (irdz.eq.0)  then 

call  ztot(nm,x,y4a,ib,p,w,zmat,ida) 
write(6,71) 

7 1  format(2x.'The  impedance  matrix  has  been  fUfed*) 

if  (iwz.eq.1)  open(unit=:23,file='zmatdat',status='unknown’) 
do  10msl,nm 
do  20  n=I^ 

if  (iwz.eq.1)  write(23,*)zmat(m,n) 

20  continue 

10  continue 

close(unit=23) 

endif 

c 

if  (irdz.eq.1)  then 

if  (iwz.eq.1)  open(unit=23,nie='zmaLdat',status='unknown') 
do  30m=l,nm 
do401sl,nm 

c  read(23,*)zmat(m4) 

40  continue 

30  continue 

write(6.72) 

72  format(2x,'The  impedance  matrix  has  been  read  in') 
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close(unit=23) 


oidif 

C 

C 

c  THIS  CORRECTS  THE  DIACX)NAL  OF  THE  IMPEDANCE  MATRIX  FOR  NON 

PEC  SURFACE 
C 
C 

call  zimped(zmaux,y4a4b.za,zb.nm4da) 
c 

write(6,73) 

73  fonnat(2x.The  impedance  diagonal  has  been  adjusted*) 
c 

c 

c  THIS  PART  ENDS  THE  PROGRAM  IF  YOU  JUST  WANT  THE  IMPEDANCE 

MATRIX 

C 

if  (ngo.eq.O)  goto  998 
c 
c 
c 

c  THIS  PART  DETERMINES  THE  BISTAHC  RADAR  ECHO  WIDTH 

c 

if  (isc.eq.l)  then 
c 

call  volt(vmaUc.y.ia.ib,nm,phip,w,ida,iwv) 
write(6.74) 

74  fofmat(2x,The  voltage  matrix  has  been  computed  for  incident  angle*) 

write(6,5)phip 

c 

caU  amp(zmat4mat,vmat4un4;heck4da,iwi) 
write(6,75) 

75  format(2x,The  current  matrix  has  been  computed*) 

C 

do  SO  i=(degstart/deg),(degstop/deg) 
phi=i*deg 

call  scat(ef,phi,x,y,ia,ib,w,nm,ida) 
sums(0.,0.) 
c 

do  60  n=l4im 
sum=sum-Kimat(n)*ef(n)) 

60  continue 

c 

es=sum 

o2d=2*pi*((cabs(es))**2) 
wavelengthsi2*3. 141592654*3.0e8/w 

o2dlog=10*logl0(o2d/wavelength) 

c 

phase=atand(imag(es)/real(es)) 

c 

if  (leal(es).ltO)  phase=phase4-180 
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if  Q>hase.|tl80)  pluse=phase<360 
wnte(26,*)phi,o2dlog 
write(27.*)phi.o2d.o!Mlog,cabs(es), phase 
wri^6,76) 

76  fonnat(2x,The  bistatk  field  has  been  detennined  for  angle*) 
writ^6,5)idii 

SO  continue 
c 

endif 

c 

c  THIS  PART  DETERMINES  THE  MONOSTATIC  RADAR  ECHO  WIDTH 
c 

if  (isc.eq.O)  then 


do  SS  i=Kdegstaii/deg),(degsu^deg) 

phi=i'*deg 

phip=phi 

call  volt(vmat.x.y,ia,ib.nin,phip.w4da,iwz) 


call  ainp(zn)at4mauvmaunni,check,ida) 
check=2 


call  scat(ef,phi,x.y4a,ib,w,nm,ida.iwi) 
sums(0.,0.) 

do  6S  nsl^im 
sum=sum-Kimat(n)*ef(n)) 
continue 


es=sum 

o2d=2*pi*((cabs(es))**2) 

o2dlog=10*logl0(o2d) 

phasesatand(imag(es)/ieal(es)) 


if  (feal(es).lLO)  phase=phase-t-180 
if  (pha».gLl80)  phase=phase-360 
write(26,*)phi,o2dlog 
write(27,*)phi,o2d,o2dlog,cabs(es),phase 
write(6.77) 

77  format(2x,'The  monostatic  field  has  been  calculated  for  angle*) 
write(6,S)phi 

55  continue 
c 

endif 

c 

c 

c 

998  stop 
end 
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c 

c  END  OF  THE  MAIN  PROGRAM 

c 

c 

^^^^^0*4i0*^^^^*ti*it*4**************************************************** 

««« 

C 

44444*0***************************************************************** 

*** 

C 

************************************************************************ 

*** 

c 

4444444*4*************************************************************** 

*** 

C 

C 

c 

c 

c 

************************************************************************ 

*** 

c 

c 

subroutine  getinput(w,isc,deg4igoJireadz,p,phip,nni,x,y,ia,ib,za,zb 
1  ,ida4wz.iwi.iwv,degstart,degstop) 
c 
c 

c  This  subroutine  inputs  the  following  data 
c  f  s  frequency  in  hertz 

c  isc  s  0,1  for  mcmostatic  or  bistatic  radar 
c  deg  =:  real  number  increment  of  degrees  to  span  over 
c  degstart,stop=start  and  stop  angles  to  be  seen  in  degrees 
c  ireadz  =  0,1  tell  whether  to  compute/read  impdance  matrix  for  PEC  body 
c  iwz4wi4wv=0,l  tell  wheter  to  not  writeAvrite  imp,  volt,  amp  matrixes 
c  p  =  EV^  number  of  divisions  to  be  used  in  Simpson's  rule  integration 
c  phip  =  incident  plane  wave  angle  in  degrees 
c  np  =  number  of  node  points 

c  x(i),y(i)  s  rectangular  coordinates  of  node  points  (meters) 
c  nm  =  number  of  modes 

c  ia(i)4b(i)  s  integer  values  of  midpoint  node  points  for  each  mode 

c  za(i),zb(i)  s  complex  impdance  at  each  node  point  (ohms) 
c 
c 

integer  isc4ip,nm4a(ida),ib(ida),ireadz4igo,p 
real  wXx(ida),y(ida),deg,degstart,degstop 
complex  za(ida),zb(ida) 
c 

open(units2 1  ,file='input.dat',status='old*) 

iead(21,*)f 

w=2*3.141592654*f 

read(21,*)isc 
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o  o 


read(21,*)deg 

read(21,*)degstart 

read(21,*)<legstt>p 

read(21,*)ngo 

read(21,*)iieadz 

read(21,*)iwz 

read(21,*)iwi 

icad(21,*)iwv 

icad(21  *)p 

read(21,*)phip 

read(2iy)np 

if  (np.gtida)  then 

writ^*,*)11iefe  are  too  many  points,  reset  ida' 
writc(*,*)'ida  =  *4da 
write(*,*ynp  =  '4ip 
stop 
end  if 
c 

do  100  i=l4ip 
read(21.*)x(i).y(i) 

100  continue 
c 

read(21,*)nm 

open(unit=22,nie='graph.dat',status='unknown') 
do  110i-l4un 
read(21.*)ia(i).ib(i) 
c 

c  This  part  writes  to  the  graph  file  for  a  gnuplot  of  our  gemnetry 
c 

write(22,*)x(ia(i)),y(ia(i)) 
write(22,*)x(ib(i)),y(ib(i)) 
if  (i.eq.nm)  wtite(22,*)x(ib(i)).y(ib(i)) 
c 

110  continue 

close  (unit=22) 
c 

do  120isl,nm 
read(21,*)za(i),zb(i) 

120  continue 
c 

close  (unit?=21) 

return 

end 


C 

C 

subroutine  ztot(nm4i,y,ia,ib,p,w,zmat,ida) 
c 
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c 

c  This  subroutine  determines  the  PEC  impedance  matrix  using  ZSELF  and  2[Mirr 
c  uo  =  is  the  pormability  of  free  space 

c  dn.dm  =  length  of  each  mode  (must  be  less  than .  1  wavelengths) 
c 

c 

integer  ia(ida)4b(ida).p,nm 
real  uo.w.x(ida),y(i^)Ji,dn.dm 
complex  zmat(i^4da) 
c 

uo=1.256637061e-6 
do210m=l,nm 
do  200  n=l  jun 

dn=sqrt((x(ia(n))-x(ib(n)))**2-Ky(ia(n))-y(ib(n)))**2) 

dm=sqrt((x(ia(m))-x(ib(m)))**2+^(ia(m))-y(ib(m)))**2) 

h=dn/p 

c 

if  (n.eq.m)  call  zself(zmat.m,w,p,h,ida) 

if  (n.ne.m)  call  zmut(zmat,x,ym.n,p,dn,dm,ia,ib,w.ida) 

200  continue 
210  continue 
c 

return 

end 

c 

c 

*** 


cc 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


subroutine  zself(zmat,m,w,p,h,ida) 


This  subroutine  determines  the  self  impdeance  terms  of  the  PEC  matrix 
using  equivalent  endpoint  technique.  It  also  uses  10  times  the  number 
of  divisions  used  in  the  mutual  impedance  subroutine  for  better 
accuracy.  It  needs  the  subroutine  HANK  to  run. 

feo  =  is  the  equivalent  endpoint 


integer  p 

real  k,h,si,xx,uo,c 

complex  sum/eo4i0 1 4i2,hextra,zmat(ida,ida) 


p=p*10 


h=h/10 

uo=1.256637061e-6 


k=w/(3e8) 

c=(w*uo)/(6*p) 

sum=(0.,0.) 

xx=k*h 
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call  hank(xxjt0,liextra,0) 
xx=2*k*h 

call  haiik(xxJil,liextra.O) 
feo=(4.656*hO)-(3.656*hl) 
c 

sum=feo*c 

do250i=2.((iV2)+l) 

if(i.iq.((p/2)+l))si=l 

xxKi-l)*k*h 

call  haiik(xxji2,hextra,0) 

sum=sum+(c^si*h2) 

250  continue 

^  h=h*10 

p=p/10 

zmat(ni4n)=sum 

return 

end 

c 


subroutine  zmut(zmat^,y,m,n,p,dn,dm,ia,ib,w ,ida) 

c 

Q 

c  This  subroutine  determines  the  mutual  impedance  terms  of  the  PEC  matrix 
c  pran  -  distance  between  center  points  of  the  modes 
c  pmi  =s  distance  from  the  ith  portion  of  the  nth  mode  to  the  mth  mode 

c  This  subroutine  needs  the  subroutine  HANK  to  run 

c 

c 

integer  ia(ida)4b(ida),p  .  ,  ,  , 

real  x(ida),y(ida),pmn,pmi,dn,dm,deltx,delty,si,c,w,uo,xm,ymJc,xx 
complex  sum  Ji2  Jiextra^at(ida4da) 

k=^/(3e8) 

uo=1.256637061e-6 

c=(w*uoy(12*p) 

sums(0.,0.) 

xm=(x(ia(m))+x(ib(m)))/2 

ym=^(ia(m)>fy(ib(m))y2 

xn=(x(ia(n))+x(ib(n)))^ 

yn=^(ia(n))+y(ib(n))y2 

pmn=sqrt((xm-xn)**2+(yni-yn)**2) 

c 

if  (pmn,lL(100*dn))  goto  230 
if  (pmn.lL(100*dm))  goto  230 
xxsk^pmn 
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call  haiik(xxji2,hextra,0) 
zmat(in,n)=((w*uo)/4)*h2 
goto  240 
c 

230  do220i=l,(p+l) 
si=3.+(-l)**i 
if  (i.eq.l)  si=l. 
if  (i.eq.(^l))  sial. 
deltxKx(ib(n))-x(ia(n))Vp 
delty=^ib(n)>y(ia(n))yp 
xi=x(ia(n)>Ki-  l)*deltx 
yi=:y(ia(n))+(i-l)*delty 
pnu=sqrt((xni-xi)*  *2+(ym-  yi)**2) 
xx=k*pmi 

call  ha^(xx  Ji2,hextra,0} 
sum=sum+(si*h2) 

220  continue 

zmat(m,n)=c*sum 

c 

240  return 
end 
c 
c 

* 

C 

C 

subroutine  zimped(zmaui.y4a,ib,2a,zbjun,ida) 
c 
c 

c  This  subroutine  determines  the  diagonal  impdeance  matrix  terms  used  to 
c  conect  our  impdeance  matrix  for  non  PEC  surfaces.  Assuming  piecewise 
c  linear  impedances, 
c 
c 

integer  ia(ida)4b(ida)4un 
real  x(ida),yOda),dn 

complex  zmat(idtuida),za(ida),zb(ida),zmn  j 

j=(0.,l.) 
do  50  i=l,nm 
zinn=(za(i)+zb(i)V2 

dn=sqtt((x(ia(i)yx(ib(i)))**2+(y(ia(i))-y(ft^^^^ 
znm=zmn/dn 
zmat(i4)=zmat(U)+zmn 
SO  continue 
c 

return 

end 

c 

4»l>* 
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c 

c 

c 

c 

*** 

subroutine  volt(vmaU(,y4a4b4un,pltip,w,idaJiwv) 
c 
c 

c  This  subroutine  determines  the  voltage  matrix  from  a  plane  wave, 
c  phip  =  incident  plane  wave  angle  in  Egress 
c 

integer  ia(ida)4b(ida)4mn.c 
real  xl,yl,phip.w,x(ida),y(ida) 
complex  vmat(ida)  j 

j=(0..1.) 

c=3e8 

c 

if  (iwv.eq.l)  open(unit=24,file=*vmat.dat‘^tatus='unknown’) 

do  300  l=l,tun 

xl=(x(iaa))+x(iba))y2 

yl=(y(iaa))+y(ib0))y2 

vmat0)=«xp(j*(w/c)*((xl*cosd(phip))4<yl*sind(phip)))) 
if  (iwv.eq.l)  write(24,*)l,vmat0) 

300  continue 

close  (unita24) 
return 
end 
c 
c 
c 

*** 

c 

subroutine  amp(zmat,imat,vmat,nm,check4da4wi) 
c 
c 

c  This  sutxoutine  determines  the  current  matrix  from  die  impedance  and 
c  voltage  matrices.  It  also  destroys  the  impedance  matrix  in  the  process 
c  This  sul»outine  needs  the  CROUT  subroutine 
c  imat = currait  matrix 
c  iwi=0,l  to  not  or  write  the  current  matrix 
c 

integer  nm, check 

complex  zmat(ida.ida),vmat(ida)4mat(ida) 
c 

if  (iwi.eq.l)  open(unit=25,file='imat.dat',status='unknown') 

do  43  l=l4im 

imat(l)=:vmatG) 

43  continue 

call  crout(zmat,imat,ida,l,0,check4un) 
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do  45 Isl^un 

if  (iwi.eq.l)  write(25,*)l4niat(l) 

45  continue 
c 

close  (units25) 
return 
end 
c 
c 

*** 

subroutine  scat(ef,phi,x,y,ia.ib.wj)ni4da) 
c 
c 

c  This  subroutine  determines  the  scattered  field  from  a  point 
c 

integer  ia(ida)4b(ida) 

real  phi.x(ida),y(ida),alpha.dn,k,deltx,delty,uo,pi,psi,a 
complex  efCida)  j  Jambd^c 

^  uo=1.256637061e-6 

pi=3, 141592654 
lambda=:(0.,0.) 
j=(0.,l.) 
k=w/(3^) 

css-((w*uoy4)*sqrt(2/(pi*k))*exp(j*(pi/4)) 

c 

do  400  nsl,nm 
deltxs:x(ib(n))-x(ia(n)) 
deltys=y(ib(n))-y(ia(n)) 
dn=sqrt(deltx**2+delty**2) 
c 

if  (deltx.eq.O)  then 

if  (delty.gtO)  al(riia=90 
if  (delty.lt0)  alpha=:270 

endif 

c 

if  (deltx.ne.O)  then 

alpha=atand(delty/deltx) 
if  (deltxJtO)  alpha^pha+180. 

endif 

c 

c 

a=cosd(phi-alpha) 

if  (a.ne.O)  lambda=((exp(j*k*a*dn))-iy(j*k*a) 
c 

if  (a.eq.O)  lambda=dn 
c 

psi=4*((x(ia(n))*cosd(phi))+(y(ia(n))*sind(phi))) 

efOiMt^<ln)*cxp(j*psi)'*lambda 

c 

400  continue 
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u  u 


return 

end 


*** 


SUBROUTINE  HANK(XXJIJI1.ID) 

C 

C  HANKEL  FUNCTION  OF  THE  SECOND  KIND 

C 

C  XX*  REAL  ARGUMENT 
C  H  =  HOFORDERO 
C  Hl*HOFORDERl 
C  ID  *0  IMPLIES  COMPUTE  H  ONLY 
C  1  IMPLIES  COMPUTE  HI  ONLY 

C  2  IMPLIES  COMPUTE  BOTH 

C 

COMPLEX  H4H 

DOUBLE  PRECISION  D,DX.C.EAJEB.FAI^ 

X=ABS(XX) 

DX=X 

IF(X.GE.3.XiO  TO  10 
C=DX*DX/9. 

D=KDLOG(DX)-.693 147 18056)*.63661977237 
IF(ID.EQ.1)GOT08 

B  =(((((.21E-3*C-.39444E-2)*C+.444479B-1)»C-.3163866)*C+.1265 
16208E+1)*C-2.2499997)*C+1 . 

Y  *D*B  +.36746691+(((((-.24846E-3»C+.427916E-2)*C-.4261214E-l) 

2*C+.25300117)*C-.74350384)*C+.60559366)*C 

H=CMPLX(B,-Y) 

IF(ID.EQ.O)RETURN 

8  B 1  =((((((.  1 109E-4*C-.3 176 lE-3)*C+.4433 19E-2)*C-.3954289E-  1)*C+ 
3.21093573)*C-.56249985)*C+.5)*DX 

Y1  =D*B1  +((((((.27873E-2*C-,400976E-1)*C+.3123951)*C-1.3164827) 

4*C+2.1682709)*C+.2212091)*C-.6366197724yDX 

H1=CMPLX(B1.-Y1) 

RETURN 
10  D=3.0D0/DX 
C=iyDSQRT(DX) 

IF(ID.EQ.l)GOT018 

EA=C*((((((. 14476E-3*D-.72805E-3)*D+.  137237E-2)*D-.95 12E-4)*D- 

5.552740B-2)*D-.77E-6)*D+.797884560803) 

FA=(((((.13558E-3*D-.29333E-3)*D-.54125E-3)*D+.262573E-2)*D- 

6.3954E-4)*D-.4166397E-1)*D-.785398163397+DX 

B=EA*DCOS(FA) 

Y=EA*DSIN(FA) 

H=CMPLX(B,-Y) 

IF(ID.EQ.O)RETURN 

18  EB=C*((((((-.20033E-3*D+.l  13653E-2)*D-.24951 1E-2)*D+.17105E-3)*D+ 
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7.l659667E-l)*D+.156E-5)*D+.797884560803) 

FB*(((((-.29166E-3*D+.79824B-3)*D+.74348E-3)*D-.637879E-2)*D+ 
8.5650E-4)*D+. 12499612)*D-2.356194490192+DX 
BlsEB*DCOS(FB) 

Y1*EB*DSIN(FB) 

HlaCMPLX(Bl.-Yl) 

RETURN 

END 


c 

c 

c 

c 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 


SUBROUTINE  CROUT(C.S.ICX:4SYM,IWR,I12,N) 

GROUT  SOLVES  A  SYSTEM  OF  SIMUL.  EQS. 

INPUTS: 

CaJ)*COEF.  MATRIX 
SO)  =  RIGHT  HAND  SIDE  VECTCm 
ICC  =  DIMENSION  OFC  ARRAY 
ISYM  =  0  IMPLIES  Ca J)  IS  SYMMETRIC 
=  1  IMPLIES  CaJ)  IS  NOT  SYMMETRIC 
IWR  =  1 ,0  TO  WRITE  OR  NOT  WRITE  SOLUTION  VECTOR 
112  =  1  ON  FIRST  CALL  TO  CROUT  WITH  A  GIVEN  CGJ) 
=  2  ON  SUBSEQUENT  CALLS 
N  *  NUMBER  OF  EQS.  AND  UNKNOWNS. 

OUTPUT; 

SO)  =  SOLUTION  VECTOR 
C(IJ)  AN  AUXILARY  MATRIX 

COMPLEX  CaCC4CC).S(l) 

COMPLEX  F.P,SS,T 
2  FORMAT(1X,1I5,1F10.3.1F15.7,1F10.0) 

5  FORMAT(1HO) 

IFai2.NE.lX30  TO  22 
IF(N.EQ.l)S(l)=S(iyC(l,l) 

IF(N.EQ.1XjOTO  100 
IFaSYM.NE.0X3O  TO  8 
DO  6 1*10^ 

D06J=IJ^ 

6  C(J.I)=CaJ) 

8  F=C(1,1) 

DO10U=2^ 

10  C(14->=C(1,L)/F 
DO20L^2J4 
LL1>L-1 
D020I=UN 

F=caj-) 
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DOllK^lOXL 

11  F=FCa.K)*C(KJ-) 

Ca4-)=F 

IF(LEQ.I)CX)T020 

PaC(U-) 

IF(ISYM£Q.OXX)  TO  15 
F=C(LJ) 

DO  12  K:=l  JiX 

12  F=F-C(L,K)*C(K.D 
C(L4)=F/P 

GO  TO  20 
15  F=cax.) 

C(L,I)=F/P 
20  CONTINUE 
22  DO30L=lJ^ 

P=C(L,L) 

T=S(L) 

IF(L.EQ.l)GOTO30 
LLL=L-1 
D025K=1JLLL 
25  T=T^(L.K)*S(K) 

30  S(L)=T/P 
D038L=2J4 
I-N-L+1 

n=i+i 

T=S(I) 

D035Ksn^ 

35  T«T-Ca.K)*S(K) 

38  S(D*T 

IF(IWR.LE.O)  GO  TO  100 
WR1TE(64) 

CNORs.O 
DO  40 1*1^ 

SA=CABS(S(I)) 

40  IF(SA.GT.CNOR)CNOR=SA 
IF(CNOR.LE.O.)CNOR=1. 

D044Ib1^ 

SS::^(I) 

SA^ABS(SS) 

SNOR=SA/CNOR 

PH=.0 

IF(SA.GT.0.)PH=57,29578*ATAN2(AIMAG(SS)JUEAL(SS)) 
44  WRrrE(6.2)l,SNOR,SA,PH 
WRrrE(64) 

100  RETURN 
END 
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Total  Ftold  (dB)  8c«M«od  FWd  (dB) 


Scattered  Field  vs  Observation  Angle  (Capacitive  Design  *1) 


Total  Field  vs  Observation  Angle  (Inductive  Design  #1) 
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Total  FMd  (dB)  Total  FMd  (dB) 


Total  FMd  vt  Azimuth  Anglo  (Inductivo  Design  #1) 


Total  Field  vs  Observation  Angle  (Inductive  Design  #1) 
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FMd(dB/lm) 


APPENDIX  E 


FDTD  RESULTS 

This  appendix  contains  all  the  FDTD  results  on  the  modeling  of  the  second 
cj^)acitive  taper.  The  results  are  shown  as  scattered  field  vs  frequency  for  different 
angles. 

FDTD  Scattered  Field  vs  Frequency 
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FDTO  Scattered  FMd  vs  Frequency 
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FIttId  (dB/m)  Scatteiwl  FMd  (dB/m) 


FDTD  Scattered  Field  vs  Frequency 
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FDTD  Scattered  Reid  vs  Frequency 
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FDTD  Scattered  Field  vs  Frequency 
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